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Non-identical particle femtoscopy probes not only the size of the emitting system, but also 
the emission asymmetries between particles of different mass, which are intimately related 
with the collective behavior of matter. We apply the technique to the simulations from the 
THERMINATOR+Lhyquid model of the heavy-ion collisions at ^Snn = 200 GeV. We present pre- 
dictions for all pairwise combinations of pions, kaons, and protons, and discuss their interpretation. 
We show that kaon and proton distributions are strongly influenced by flow: the source gets smaller 
and shifted to the outside with growing pr, while for pions the shift is significantly smaller, produc- 
ing an emission asymmetry. We explain how particles coming from decays of hadronic resonances 
enhance the asymmetry signal coming from flow, contrary to naive expectations. Emphasis is put 
on extracting this unique information about collective behavior of matter from the non-identical 
particle correlations. We also present, in detail, the technical aspects of the non-identical particle 
femtoscopy technique applied to data from the heavy-ion collisions. We list the sources of system- 
atic errors coming from the method itself and the usual assumptions. We describe robust analysis 
methods and discuss their limitations. 
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I. INTRODUCTION 

Femtoscopy has been used for more than 35 years [l|, [1] 
to measure sizes of the systems created in nucleus-nucleus 
collisions. Initially it was developed for analyzing the 
two-particle correlation arising from the wave-function 
symmetrization for pairs of identical particles Q and was 
similar in mathematical framework to the "HBT interfer- 
ometry" used in astronomy Later it was realized 

that similar correlations arise due to the Final State In- 
teractions (Coulomb and Strong) between particles that 
are not necessarily identical |7Hlll|. 

Femtoscopy of non-identical particles provides unique 
information: while the identical particle correlation usu- 
ally only measure the "size" of the emitting region (more 
precisely, the second moments of the emission function), 
non-identical correlations can, thanks to the very fact 
that they correlate particles that are not identical, also 
measure the relative emission shifts (the first moments of 
the emission function) 0^,[l3-[ll|. 

The unique features of the non-identical particle fem- 
toscopy was used in low-energy nuclear collisions to study 
the time ordering of the emission of various nuclear frag- 
ments from the compound nucleus [T5l - [l9| . This required 
measuring time differences from several to hundreds of 

fm/c 

In this work we focus on the applications of this tech- 
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nique in the collisions of ultra-relativistic heavy ions, 
specifically Au ions at the Relativistic Heavy Ion Collider 
(RHIC) in Brookhaven National Laboratory (BNL), al- 
though the discussion is also relevant for lower (at Super 
Proton Synchrotron SPS of CERN) and higher (Large 
Hadron Collider LHC at CERN) energy collisions. 

One of the major discoveries at RHIC has been the 
observation of the Quark Gluon Plasma (QGP), which 
was found to behave collectively, very much like a 
fluid [2^ [2l|. Hydrodynamic equations seemed to de- 
scribe this behavior well in the momentum sector. How- 
ever attempts to simultaneously describe the space-time 
behavior measured by femtoscopy have not been success- 
ful until recently. This failure was commonly referred 
to as the "RHIC hydro HBT puzzle" ilH. Through a 
detailed analysis of the experimental data and various im- 
provements of the hydrodynamic description it was real- 
ized that certain assumptions had to be modified in order 
to properly describe both sets of observables. The initial 
condition was changed to the one that uses Gaussian ini- 
tial profile for the transverse energy density 26|, which 
results in faster development of the initial flow, compared 
the traditional one [13, [13 • The equation-of-state used 
did not exhibit a first-order phase transition, but rather 
a cross-over. The detailed simulation of the resonances 
contribution in the later stage of the collisions was carried 
out [13 . In addition, some studies suggest that the intro- 
duction of viscosity and universal (pre-equilibrium flow) 
into the model may also play an important role [H, [s^l • 

The hydrodynamic scenario produces specific space- 
momentum correlation patterns, which are commonly 
referred to as flow. The system created in the heavy- 
ion collision expands rapidly outwards, showing a very 
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strong radial flow, which is observed in the modification 
of the single particle inclusive px spectra shape. In ad- 
dition, in non-central collisions, the initial overlap region 
has an elongated (usually described as "almond" ) shape, 
and this spatial asymmetry is converted to the momen- 
tum one in the final state. This momentum asymme- 
try is observed as an elliptic flow V2 and is the subject 
of v ery intensive theoretical and experimental studies; 
see [sj, nil for recent reviews. However, because of its 
origin, it is small in central collisions. This is unfortu- 
nate, since it is in central collisions that we expect to 
create the largest volume of deconflned matter, which we 
would like to study. Both the pr spectra and the ellip- 
tic flow are observables depending only on the momenta 
of the particles, so their connection to space-time can 
only be indirectly inferred. To access it directly we em- 
ploy femtoscopic techniques. It is argued that the fall 
of the "femtoscopic radii" with particle's rriT can be in- 
terpreted as the decrease of "lengths of homogeneity" a 
direct consequence of radial and longitudinal flow |33| . 
However, one might come up with alternative explana- 
tions, involving temperature gradients to produce sim- 
ilar dependencies i34]- In this work we will show how 
the collective flow present in hydrodynamic and trans- 
port models, in addition to the effects discussed above, 
produces differences in average emission points between 
particles of different masses. We will also describe how 
these emission asymmetries can be accessed via the non- 
identical particle femtoscopy js^l . Measuring such effects 
would enable us to eliminate scenarios alternative to hy- 
drodynamic expansion and provide a crucial and strict 
test for the models. Results on 7r+-7r~ and 7r-proton cor- 
relations have been previously reported by NA49 [35, 36] 
and the CERES experiment ^[si] at the SPS. The first 
results on pion-kaon correlations have been reported by 
STAR [3^ . Preliminary results on tt-S correlations from 
STAR have also been presented [i^. Proton-A correla- 
tions have also been measured by STAR [4l|, but this 
analysis did not attempt to extract emission asymme- 
tries. 

In order to properly simulate the emission asymme- 
tries between particles of different masses we need a 
model which has all the important features: the hydro- 
dynamic phase which produces space-momentum correla- 
tions and the hadronic phase where at least the hadronic 
resonance decays and propagation is treated. In addi- 
tion the model should provide the space-time freeze- 
out coordinates of particles, so that femtoscopic calcu- 
lations can be carried out. It should, as much as pos- 
sible, reproduce the available data on particle spectra 
and femtoscopy. Hydrodynamics inspired "blast-wave" 
parametrization of freeze-out have been used to model 
asymmetries between non-identical particles j43 |. This 
study nicely illustrated the connection between strong 
collective matter behavior and the asymmetries; however 
it used simplified emission functions and neglected the 
resonance propagation and decay. In this work we have 
chosen the THERMINATOR+Lhyquid model, which we in- 



troduce in Section ini In Section Hill we describe the par- 
ticular set of simulations of heavy-ion collisions at six 
centralities that we have carried out. We calculated the 
system size and emission asymmetry for three pair types: 
pion-kaon, pion-proton and kaon-proton. We discuss 
the origins of the asymmetry and provide the connec- 
tion to the space-momentum correlations coming from 
fiow. We also discuss other non-flow sources of emission 
asymmetries and provide a quantitative estimate of all 
contributions. 

Having emphasized the importance of the emission 
asymmetries phenomenon we proceed to describe the 
theoretical framework of the non-identical particle cor- 
relations in Section IIVI Then, we move to the techni- 
cal aspects of the measurement in heavy-ion collisions 
in Section |Vl We argue that a specific mathematical 
representation of the non-identical correlation functions 
(the spherical harmonics decomposition) shows remark- 
able synergies with the analysis and maximizes the statis- 
tical significance of the emission asymmetry signal. We 
describe the measurement procedure which enables to re- 
cover the properties of the model emission function from 
the "experimental" correlation function. In Section IVII 
we test the robustness of the procedure, list the assump- 
tions and the approximations that need to be employed 
and estimate the systematic error coming from the tech- 
nique itself. We show that the most important contri- 
bution to this error is the correct estimation of the frac- 
tion of non femtoscopically correlated pairs (tradition- 
ally called "purity" ) . We analyze this effect in detail and 
present purity estimates, based on our model of choice, 
which can be directly used in the experimental analysis. 

Finally, in Section IVIII we use numerical simulations 
with the THERMINATQR model to produce "experimental- 
like" correlation functions, which we then analyze with 
the methods presented earlier. We demonstrate that the 
method is able to recover the theoretical input values. 
The results presented in this section are the theoreti- 
cal predictions for the overall system size as well as the 
emission asymmetry for all considered pair types and 6 
centralities. They can be immediately compared to the 
experimental results with minimal number of approxima- 
tions, and therefore minimal systematic uncertainty. 



II. LHYQUID + THERMINATOR MODEL 

In this work we use what we call the standard ap- 
proach, consisting of ideal-fluid hydrodynamics followed 
by statistical hadronization. Numerous calculations have 
been performed in this framework, with the common dif- 
flculty of simultaneously describing femtoscopy and 
other signatures in the data. More precisely, the RHIC 
HBT puzzle [22| - [25| refers to problems in reconciling the 
large value of the elliptic flow coefficient, V2, with the 
identical particle interferometry in calculations based on 
hydrodynamics (43l - |47j . Recently, a successful uniform 
description of soft observables at RHIC has been ac- 
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complished, including the femtoscopic radii, within the 
standard approach [2g|. The essential ingredients of this 
analysis are the Gaussian initial condition for hydrody- 
namics, early start of the evolution, the state-of-the art 
equation of state with smooth crossover, and the use of 
THERMINATOR ^ with all resonances from SHARE 
incorporated to carry out the statistical hadronization 
at the freeze-out surface of temperature 150 MeV. The 
interplay of these elements resulted in a simultaneous de- 
scription of the transverse-momentum spectra of pious, 
kaon and protons, the V2, and the femtoscopic correlation 
radii of pions, including full centrality, kx and reaction 
plane dependence of azimuthally sensitive HBT signa- 
tures [501 • 

In this section we describe the essential elements of our 
method to the extent they are necessary for the presented 
new results. More details concerning the hydrodynam- 
ics can be found in Refs. [1^, [5l|, [l^l, while the method 
used for femtoscopic calculations has been presented in 
Rcf. M- 



A. Initial condition 

As reported in Ref. [1^ , the use of the initial condition 
for hydrodynamics of the Gaussian form, 



n(a;,y) = exp( ^ 



y_ 

262 



(1) 



where n is the initial energy density, while x and y de- 
note the transverse coordinates, leads to a much better 
uniform description of the data for the pT-spectra, V2, 
and the pionic femtoscopic radii compared to the use of 
the standard initial condition from the Glauber model. 

The width parameters a and 6 depend on centrality. 
In order to estimate realistic values for them we run the 
GLISSANDO [Eil Glauber Monte Carlo simulations which 
include the eccentricity fluctuations [5^ [5^. Then we 
match a2 and 62 to reproduce the values (a;2) and {y'^) 
from the GLISSANDO profiles. Thus, by construction, the 
spatial RMS radii of the initial condition and its eccen- 
tricity is the same as the ones from the Glauber calcu- 
lation. Nevertheless, the shape is not the same. The 
Gaussian profiles are sharper near the origin, which re- 
sults in a faster buildup of the Hubble-like flow in the 
hydrodynamical stage. 

The Glauber calculations, needed to obtain the a 
and 6 parameters, correspond to the mixed model (57| . 
where the number of produced particles is proportional 
to (1 — a)Nyj/2 + aiVbin, with iVy, and iVbin denoting the 
number of wounded nucleons jHl] and binary collisions, 
respectively. The parameter a = 0.145 for top RHIC en- 
ergy p9, 60\. The inelastic nucleon cross section is 42 mb 
for RHIC [6l|. The simulations incorporate the fluctua- 
tions of orientation of the fireball (the variable-axes ge- 
ometry) , which result in increased eccentricity compared 
to the fixed-axes geometry [g^]- Finally, the expulsion 
distance of 0.4 fm is used in the generation of the nuclear 



distributions. A source-dispersion parameter of 0.7 fm is 
used. It describes the random displacement of the source 
from the center of the wounded nucleon or the binary- 
collision position [545. 

The values of the a and 6 parameters for various cen- 
tralities and the corresponding eccentricity parameters 



(2) 
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are collected in Table ID 

The energy-density profile ([T]) determines the initial 
temperature profile via the equation of state 52]. The 
initial central temperature, T^, is a parameter dependent 
on centrality. For RHIC calculations it is adjusted to 
reproduce the total particle multiplicity. 



B. Hydrodynamics 

The hydrodynamics equations used in this work were 
described in greater detail in Refs. fsi], [s^l- We use in- 
viscid (ideal- fiuid) , baryon-free, boost- invariant hydro- 
dynamics. The equations are written in terms of the 
velocity of sound, c^, whose temperature dependence en- 
codes the full information on the equation of state of the 
system. We incorporate the known features of Cs(T), 
which, at high terriperatures, are given by the lattice 
QCD calculations [6j|; at low T they follow from the 
hadron gas including all resonances, while at intermedi- 
ate T an interpolation is used. No sharp phase transi- 
tion, but a smooth cross-over, is built in, in accordance 
to the present knowledge of the thermodynamics of QCD 
at zero baryon chemical potential. The plot of the result- 
ing Cs (T) can be found in Ref. [2^ . 

The initial proper time of the start of hydrodynamics 
is fixed to the value 



To = 0.25 fm. 



(3) 



This early start of hydrodynamics allows for a fast gen- 
eration of transverse flow, similar to the effect described 
in M. 



C. Freeze-out 

The hydrodynamic evolution proceeds until freeze-out 
occurs, where the assumed condition for the universal 
freeze-out temperature is Tf = 150 MeV. This value is 
somewhat lower than in several fits of the chemical freeze- 
out (63 - l66| : however, it agrees with the recently made 
global fits to particle transverse momentum spectra of 
Ref. '6^, Isl], where the value around 150 MeV was ob- 
tained for the kinetic freeze-out temperature.^ 



^ The use of this lower freeze-out temperature needs the intro- 
duction of the strangeness inequilibrium factors 7s in order to 
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TABLE I: Shape parameters a and b of Eq. ([T} for various centrality classes obtained by matching (x^) and (y^) to GLISSANDD 
simulations, the variable-axes eccentricity e* , and the central temperature Ti. 



c[%] 


0-5 5-10 10-20 20-30 30-40 40-50 50-60 60-70 70-80 


a [fm] 
b [fm] 
e* 

Ti [MeV] 


RHIC Au-I-Au at ^Sivjv = 200 GeV 
2.70 2.54 2.38 2.00 1.77 1.58 1.40 1.22 1.04 
2.93 2.85 2.74 2.59 2.45 2.31 2.16 2.02 1.85 
0.08 0.12 0.18 0.25 0.31 0.36 0.41 0.46 0.52 
500 491 476 460 429 390 344 303 261 



THERMINATDR is used to carry out the statistical 
hadronization at the freeze-out hypersurface according 
to the Cooper- Frye formulation [6!||. According to the 
assumed single-freeze-out approximation, identifying the 
kinetic and chemical freeze-out temperatures, rescatter- 
ing processes after freeze-out are neglected. We have 
checked that the collision rate after freeze-out is mod- 
erate for the hypersurfaces applied in this work. We es- 
timate it by considering a pion straight-line trajectory 
and counting the number of encounters with other parti- 
cles closer than the distance corresponding to the pion- 
hadron cross section. The average number of these tra- 
jectory crossings is about 1.5-1.7 per pion. This shows 
that the single-freeze-out approximation [65| works rea- 
sonably well for the present case. At a more detailed 
level, one could use hadronic afterburners to model the 
elastic collisions [70- , or attempt the hydro- kinetic ap- 
proach implemented in [23, [TjI ■ 



D. Two-particle femtoscopy 

The method used for femtoscopic analysis of the 
THERMINATOR model output was described in detail in 
Ref. 53] . The features of the formalism specific to the 
non-identical particle correlations are described in Sec- 
tion IIVI Their discussion is one of the main points of 
this work. 



III. MODEL PREDICTIONS FOR RHIC 

In this section we discuss the general trends that are 
expected to emerge in the non-identical correlations. We 
also discuss their physical origins and the possible physics 
conclusions that can be drawn when they are observed 
experimentally. Later in this paper, we will see how these 
trends manifest themselves in the calculations and we 
test whether it is indeed possible to observe them in the 
experiment. 

Following convention, we assign symbol x for the emis- 
sion points of single particles, usually used as a vector. 



reproduce the abundances of strange particles [67ll . 




FIG. 1: Diagram of the particle's velocity j3 decomposition 
into the flow /?/ and thermal f5t components. 



X, or as a vector magnitude, x. Momentum of a particle 
is denoted as p. Single particle velocity is /3. We also 
define the pair variables. Relative separation between 
particles (vector) is: r ^ xi — X2. We use the out-side- 
long coordinate system, where the long or longitudinal 
direction is along the beam axis, the out or outward di- 
rection is along the pair total transverse momentum and 
the side or sideward is perpendicular to the other two. 
In Longitudinally Co-Moving System (LCMS) the pair 
longitudinal momentum vanishes: pi^iong = ~P2,iong- 

In 

the Pair Rest Frame (FRF) (also called Pair Center of 
Mass (PCOM)) the pair center-of-mass rests: pi = ~P2- 
In our convention all pair variables in PRF are marked 
by an asterisk. The pair relative momentum half is de- 
noted as k* and is equal to the first particle's momentum 
in PRF. The pair total momentum is denoted as P (or 
sometimes K traditionally used in identical particle fem- 
toscopy). All single-particle coordinates as well as pair 
relative variables use lower case letters. In contrast we 
will use upper-case letters to denote the parameters of 
distributions. In particular we call the variance or the 
two-particle separation distribution R. We call the mean 
of such distribution fi (with the L superscript if they are 
defined in LCMS). We refer the reader to Appendix El for 
explicit mathematical formulas as well as relations be- 
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FIG. 2: Emission points of particles with the velocity in 
(0.6,0.8), its direction pointing "upwards", for kaons, central 
Au-|-Au. Panels a), c), e) show standard simulation results, 
panels b), d), f): with flow correlation "scrambled" (see text 
for details). Upper plots a), b) show primordial particles, cen- 
ter c), d): particles coming from resonance decays, lower e), 
f): all. 



tween source characteristics in the two reference frames. 

The discussion below is based on an example calcu- 
lation from the THERMINATDR+Lhyquid model, done for 
parameters tuned to the central (0-5%) RHIC Au-|-Au 
collisions at ^^snn — 200 GeV. Whenever we mention 
RHIC Au-f Au calculations we mean simulations at the 
top RHIC energy. 



A. Emission asymmetries 

Hydrodynamic evolution of matter implies strong 
space-momentum correlations in particle emission. Par- 
ticles emitted from a given fluid cell will have a velocity 
which is a combination of two components: the fluid cell 
velocity (3 / (taken from the flow field u^{r)) and the ther- 
mal velocity /3t . This is schematically shown in Fig. [TJ 

The system created in heavy-ion collision, when med- 




io 

FIG. 3: Emission points of particles with the velocity magni- 
tude in the range (0.6,0.8), its direction pointing "upwards", 
for pions, central Au+Au. Panels notation the same as in 
Fig-H 



eled in hydrodynamics, naturally develops a collective 
behavior - the radial flow in the transverse plane, that is 
matter is collectively moving "outwards" from the cen- 
tral axis of the source to the outside. In essence this is an 
X — p correlation - the direction (j)f oi the fluid element's 
transverse velocity is aligned with it's transverse position 
vector direction (j)^. When the fluid element emits parti- 
cles, all of them will have the same common flow velocity 
/3f taken directly from u^. To this velocity one adds 
a thermal component /3t which has a random direction 
(j)t in the rest frame of the fluid element. It will dilute 
the x — p directional correlation. Let us now consider the 
mean emission point of a single particle, more specifically 
its component parallel to the velocity j74j: 



^ r{l3f+l3tCos{(l>t-(f>f)) 

(3 p 



(4) 



We analyze its average over particles at fixed /3. If we 
assume a Gaussian density profile with radius Tq and 
linear transverse velocity profile /3/ = Por/rQ then we 
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TABLE II: Single-particle source parameters (mean (r) and RMS X) for three particle types in selected kinematic regions. 



Particle 


/3 


PT [GeV/c] 


Xout [fm] { 


Xout) [fm] 


Xside [fm] 


(xside) [fm] 


primordial 


K 


0.6-0.8 


0.42-0.56 


2.68 


4.46 


3.25 


0.03 


TT 


0.6-0.8 


0.12-0.16 


3.47 


2.85 


3.66 


-0.02 




0.95-0.97 


0.42-0.56 


2.56 


4.81 


3.31 


-0.01 


P 


0.6-0.8 


0.80-1.06 


2.02 


5.40 


2.82 


-0.05 


TT 


0.6-0.8 


0.12-0.16 


3.47 


2.85 


3.66 


-0.02 


TT 


0.985-0.991 


0.80-1.06 


1.88 


5.68 


2.90 


0.00 


all 


K 


0.6-0.8 


0.42-0.56 


2.99 


4.97 


3.55 


0.00 


TT 


0.6-0.8 


0.12-0.16 


4.52 


1.73 


4.72 


-0.03 


TT 


0.95-0.97 


0.42-0.56 


2.94 


5.03 


3.20 


-0.01 


P 


0.6-0.8 


0.80-1.06 


2.31 


6.12 


2.88 


-0.02 


TT 


0.6-0.8 


0.12-0.16 


4.52 


1.73 


4.52 


-0.02 


TT 


0.985-0.991 


0.80-1.06 


2.19 


6.20 


3.08 


-0.01 



TABLE III: Mean emission points and pair asymmetries for pions, kaons and protons with velocity (0.6, 0.8) in central Au+Au 
simulation. Averaging is done over all pairs in this range. See text for explanation of "flow" and "scrambled". 







{-^OTit ) 


i^out) 






\-^out 




flow 


primordial 


2.83 


4.47 


5.61 


-1.64 


-2.78 




-1.14 


non-primordial 


1.34 


7.35 


9.19 










all 


2.00 


5.54 


6.69 


-3.54 


-4.69 




-1.15 


scrambled 


primordial 


-0.04 


0.00 


-0.03 


-0.04 


-0.01 




0.03 


non-primordial 


0.88 


3.17 


4.20 










all 


0.48 


1.20 


1.28 


-0.72 


-0.80 




-0.08 



obtain \74 



/3g + T/mt ' 



(5) 



where we have exphcitly given the formula for the ve- 
locity component coming from temperature. Pions and 
kaons emitted from the fluid element will have the same 
common flow velocity /3/ (and different momenta). They 
will also get a random thermal kick, depending on mo- 
mentum. For a pion the same px kick will mean much 
larger velocity j3t kick, than for a kaon. The final ve- 
locity direction of a pion will be, on average, less cor- 
related with its emission position than that of a kaon. 
We assume that the spatial characteristics of pion and 
kaon emission are the same and the flow velocity is 
also the same. Hence, for both pions and kaons, (r/3/) 
is the same. The only difference is the T/ttit com- 
ponent, smaller for kaons. Therefore, (xout) for pions 
is smaller than {xoutj for kaons. Summarizing: when 
correlating pions and kaons with the same velocity (or 
more generally two non-identical particles of different 
masses, but same velocity) pions (lighter particles) ap- 
pear to be, on average, emitted closer to the center of 



the system than kaons (heavier particles) . Consequently, 
hydrodynamics predicts a negative emission asymmetry 

liqht. heavy / linht,heavy\ / linht heavyX i , 

' - ' rjf " ) = ( xJ,t - x„„, " ) between 



f^out ~ out I \^out -^out 

non-identical particles of different masses^. 

Let us discuss various limits in Eq. ([5]) . When there is 
no flow (/3f =0), all average emission points are zero and 
the asymmetry vanishes. If the temperature is very large 
compared to the flow velocity (or more generally the ran- 
dom component dominates over the correlated one) the 
average emission point is zero (center of the source). If 
this happens for both particles, the asymmetry is zero. 
If it happens for only one of them, the asymmetry ex- 
ists. If the flow velocity strongly dominates over temper- 
ature, and particles are emitted from the same system, 
both average emission points are strongly shifted by the 
same amount and in consequence the asymmetry is again 
small. From the discussion we see that the existence of 
emission asymmetry is not trivial and only arises in a sys- 



^ We have chosen the convention to always take the lighter particle 
as first in the pair. 
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FIG. 4: Emission points of particles with the velocity in 
(0.6,0.8), its direction pointing "upwards", for protons, cen- 
tral Au+Au. Panels notation the same as in Fig. [21 



tern where both random(thermal) and correlated(flow) 
velocities exist and are comparable in magnitude. 

We will illustrate the consequences that this mecha- 
nism has for the particle's emitting regions using the 
THERMINATOR calculations for central Au-|-Au colhsions 
at the top RHIC energy as an example. First, we focus 
on the "primordial" particles, that is particles coming 
directly from the hydrodynamical stage. Let us consider 
two variables: 



Xout = xpr/pT 

Xside = X X Pt/PT- 



(6) 



The first one is the component of the particle's trans- 
verse emission vector parallel to the particle's transverse 
momentum; the other is the perpendicular component. 
In the upper-left panel of Fig. [5] the distribution of these 
emission components are shown for primordial kaons. A 
strong correlation is seen - all particles moving "upwards" 
(to the positive "out" direction) are emitted from the 
positive Xout part of the source. This focusing of par- 
ticle's emission has two distinct effects: the overall size 



of the emitting region shrinks both in the "out" and the 
"side" direction and the average emission position shifts 
in the "out" direction, but not in "side" . Looking at the 
upper-left panel of Fig.[3]one sees the same effects but to 
a much smaller degree for pions with the same velocity: 
the size is shrunk, but not as much, the average emission 
position is also shifted but by a smaller amount. The 
exact numerical values for the "size" (or variance) and 
"shift" (or mean) of these distributions can be found in 
Tab, im Fig. 13] shows analogous pictures for protons. As 
expected, the effect is the strongest here, as protons have 
almost twice the mass of kaons. 

From the discussion above, and the numerical values in 
Tab.lUand lllll one immediately sees that hydrodynamics 
produces two distinct trends in femtoscopic observables: 
(a) The size of the emitting system gets smaller with 
increasing px of the particle. This effect is well under- 
stood theoretically (so called "lengths of homogeneity" ) 
and universally observed experimentally in femtoscopy 
in heavy-ion collisions ("tot scaling") [2^. (b) Average 
emission points of particles with different px (for exam- 
ple with same velocity but different mass) are different, 
and this difference is well defined: lighter particles ap- 
pear to be, on average, emitted closer to the center of 
the source"^. 

The second trend is the main focus of non-identical 
particle analysis. It so happens, that the non-identical 
particle femtoscopy correlates particles with the same ve- 
locity but possibly with different masses. Moreover, it 
has a unique feature of being able to measure not only the 
"size" (more precisely - the second moment of the two- 
particle distribution - the variance) but also the "shift" 
(i.e. the first moment of the distribution - the mean) 
between average emission points. It is therefore able to 
directly test the predictions of the hydrodynamic model 
with respect to the x — p correlations. We would like 
to emphasize that this is the most direct and unambigu- 
ous signal of collectivity available to femtoscopy. The tot 
scaling, which is predicted by hydrodynamics, can also be 
explained by other mechanisms not requiring collectivity 
(e.g. "temperature gradients" [s^l); but no mechanism 
is known, which would produce such specific emission 
asymmetries with no collectivity. 

The asymmetry is predicted to arise only in the "out" 
direction. The "side" asymmetry is zero. It can also be 
shown that for rapidity symmetric systems of a collider, 
such as RHIC, when the target is identical to the pro- 
jectile, the longitudinal asymmetry is also expected to 
vanish. 



^ Note that the values in Tab [ll] are calculated vs. the single- 
particle momentum direction, which is only the approximation 
of the out direction. Values in Tab. IIIII are calculated vs. the 
proper out direction - the total momentum of the pair. Therefore 
small differences between the shift values in the two tables are 
to be expected. 
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B. "No directional correlation" test 

In the last section we have argued that the correlation 
between spatial emission angle (f)r and particle's velocity 
direction (j)f is responsible for the emission asymmetries. 
We have tested this argument by performing a calculation 
in which we have intentionally broken this correlation, in 
order to show that in such case no asymmetries arise. We 
take each primordial particle separately. From its original 
transverse emission coordinates {xo,yo) we calculate its 
transverse emission radius r and angle 4>r- Then, we 
randomize the angle (pr, but keep the emission radius r 
unchanged, and calculate the new "scrambled" emission 
point {xs^ Us)- If the particle is unstable and consequently 
decays, the emission points of all daughter particles are 
shifted by the same amount, calculated for the parent 
particle (xs — Xo,ys — Ho)- Note that the momentum 
observables are not affected. 

The upper-right panels of Figs. [2 [3] and |4] show the 
effect of this procedure on primordial particles. As ex- 
pected, the average emission point of all particle types 
is now at (0, 0), and the asymmetry between particles of 
different masses, shown in Tab. IIIII is zero. Our claims 
are confirmed: the model with no directional correlation 
shows no asymmetry. 



C. Importance of resonances 

The Lhyquid+THERMINATOR model assumes that the 
evolution of the heavy-ion collision proceeds in stages. 
After the initial non-equilibrium phase there is a collec- 
tive phase, which can be well described by hydrodynamic 
equations, with the equation of state assuming the exis- 
tence of the Quark Gluon Plasma. However, at some 
point the system becomes so dilute that the continuous 
medium description of hydrodynamics is no longer vi- 
able and one converts the system to hadronic degrees of 
freedom - the THERMINATOR model does it via the von- 
Neumann sampling of the probability distribution ob- 
tained from the Cooper-Frye formula. All known reso- 
nances are the degrees of freedom in the hadronic phase. 
Their abundances are well described by chemical models, 
and if one trusts these calculations, then at least 2/3 of 
the observed pions do not come from the original hydro- 
dynamic phase (the so-called "primordial" particles) , but 
are daughters of resonances. Additionally all resonances 
have their intrinsic lifetime, so they travel some distance 
before decaying - hence they will certainly modify the 
space-time picture of particle emission. It is clear that a 
careful and detailed simulation of the resonance propaga- 
tion and decay (e.g. as implemented in THERMINATOR) is a 
critical feature of the model which aims to describe fem- 
toscopic observables. Also, the relative abundances and 
decay momenta of resonances producing pions, kaons and 
protons as final particles are obviously quite different; 
so, the feature is even more important for non-identical 
particle correlations, where we study relative differences 



between various particle types. 

Let us first qualitatively consider a resonance decay 
process. We have the original resonance, produced as a 
primordial particle. Its emission followed common "flow" 
X — p correlation, as for any other particle. Resonances 
are usually quite heavy, so we expect this correlation to 
be strong. The resonance travels some distance with the 
original velocity - this enhances the x — p correlation. 
After some random time the decay process occurs. Two 
(or three) particles are created at the decay point. We 
note the similarity between the resonance decay process 
and the emission from the fluid element. The daugh- 
ters of the resonance will have the "common" velocity - 
this time it is not the fluid elements velocity, but simply 
the resonance's one. And they will have the "random" 
component - the "decay momentum" of the given decay 
channel. And again the random component will matter 
more for the lighter particle - the velocity that corre- 
sponds to the fixed "decay momentum" will be larger for 
the lighter particle. One can imagine two scenarios. If 
the decay momentum of the given decay channel is large, 
compared to the daughters' mass, than the daughters' 
emission direction will be randomized, and the common 
x—p correlation will be lost. On the other hand if the de- 
cay momentum is small, the correlation will be preserved 
(or, in other words, the daughter particle will travel in 
roughly the same direction as the parent) or even en- 
hanced due to the additional x—p correlation resulting 
from resonance propagation. 

From the description above one concludes that the res- 
onance process will induce the space-momentum correla- 
tion in a way similar to the collective flow; so, one can 
ask if it can be an alternative mechanism producing such 
asymmetries. However, qualitative expectations for val- 
ues of the asymmetries or even general trends are not 
immediately obvious. They will non-trivially depend on 
relative abundances of resonances and their daughters, as 
well as particular values of decay momenta in specific de- 
cay channels. The fact that some resonances decay in cas- 
cades makes it even more complicated . Moreover the 
original primordial resonances will have a natural x — p 
correlation coming from the earlier hydrodynamic phase. 
One then faces a quantitative problem: do the resonance 
decays introduce emission asymmetries on their own, in- 
dependent from flow asymmetries? And if yes, how big 
are they, compared to the asymmetries coming from flow? 
In particular are they small enough, so that one can still 
safely interpret the asymmetry observed in the experi- 
ment as coming from collective behaviour such as flow? 

To answer these questions a detailed simulation is 
needed, in particular one that has intrinsic x — p correla- 
tions coming from flow (preferably with the possibility to 
switch them off) and which incorporates all known res- 
onances, together with the state-of-the art knowledge of 
their masses and decay channels. From this description 
it is clear that THERMINATOR model is perfectly suited for 
the task. In addition, we use the "scrambling" procedure 
described in the previous Section to switch off the x — p 
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correlations coming from flow, in order to estimate the 
asymmetry coming from the decay processes alone. 

Let us see what is the effect of resonance decays on 
the properties of the emitting regions. In panels e) of 
Figs. (sunn] the emission points of all particles (both pri- 
mordial and from resonances) are shown, while Tab. HIl 
shows the numerical values of sizes and shifts. The re- 
sulting asymmetries are listed in Tab. IIIII We consider 
the size of the system first. As expected the overall size is 
larger - that is naturally expected as the resonances travel 
some distance before they decay. The mr scaling seems 
to be preserved. In general the trend is consistent with 
previous studies of resonance influence on femtoscopic 
observables [H, 1?^ and agrees with the natural expecta- 
tions. Less trivial and more interesting effects are visible 
in the average emission points. A qualitatively different 
effect is seen for pious than for kaons and protons. We 
first inspect the emission points for particles coming from 
resonances only (no primordial ones), seen in panels c) 
of Figs. [HI [21 m For pious the average emission point is 
shifted less from the center for resonance daughters, than 
for the primordial ones. Apparently the first scenario 
described three paragraphs before is in effect here: the 
decay momenta of the resonances producing pious are so 
large, compared to the pion mass, that they completely 
wash out the original fiow x — p correlation. In contrast 
both for kaons and protons the resonance daughters are 
shifted more from the center than primordial particles. 
This time the second scenario is in effect - the decay mo- 
menta are small compared to particles' masses, so small, 
in fact, that they are not even able to counter the ad- 
ditional X ^ p correlation coming from resonance prop- 
agation. A more detailed discussion of this effect, with 
examples of particular resonances and decay channels for 
pious and kaons is found in . The effect persists when 
one takes all particles, primordial and resonance daugh- 
ters, together. Both effects collaborate in enhancing the 
pion-kaon and pion-proton asymmetries: 

(^out) ~ i^out) ~ i^out) 

{'"'out) ~ i^out) ~ i^out) C^) 

while the kaon-proton asymmetry stays rather similar to 
the primordial only case (and small) : 

We compare the source distributions in panels on the left 
in Figs. [21 [H [Uto the ones in panels on the right, where 
the "scrambling" procedure was applied. The numeri- 
cal values for asymmetries are given in Tab. IIIII First, 
let us focus on pion-kaon and pion-proton pairs. As al- 
ready discussed, "scrambled" primordial particles show 
no asymmetry. The non-primordial ones do show some, 
but still significantly lower than the non-scrambled ones. 
This shows that for non-primordial particles both sources 
of asymmetry are important: the original x — p correla- 
tion of the parent particle and the additional asymmetry 



TABLE IV: Time and space asymmetries (from the fits to 
the distributions around their peaks) for central Au+Au col- 
lisions. PRF values are with asterisk. 
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from the decay process. However, the former dominates. 
A critical test is the comparison of asymmetries for all 
particles between normal "fiow" and "scrambled" scenar- 
ios, since this is the observable measured in the exper- 
iment. As the simulation shows, the asymmetry which 
can be attributed solely to the trivial resonance decay 
processes can account for only 20% (17%) of the total 
asymmetry produced in central Au-|-Au collisions at top 
RHIC energy for the pion-kaon (pion-proton) pair. For 
kaon-proton the additional asymmetry produced by the 
resonance decays is negligible. 

In summary, we have shown that resonance decays do 
not dilute, but rather enhance the asymmetry signal for 
pion-kaon and pion-proton pairs. At the same time, even 
though the resonance decay process can potentially be 
an independent source of emission asymmetry, we have 
shown, by detailed calculations, that in realistic condi- 
tions such asymmetry is less than a quarter of that pro- 
duced by flow. Therefore our original expectation holds: 
if significant emission asymmetry is observed in pion- 
kaon and pion-proton correlations in the experiment, it 
favors the explanation of strong x — p correlations in the 
emitting system, such as the ones produced by hydrody- 
namic radial fiow. 



D. Flow versus time asymmetries 

As we have noted in the Introduction, the non-identical 
particle correlations technique has been initially devel- 
oped to measure mainly the emission time differences be- 
tween various particle species. Only later the connection 
to radial flow was noted [s^ . 

In Au-t-Au collisions one expects intrinsic time differ- 
ences^ for at least two reasons. The reader is referred 
to -SQl , where the evolution of the average emission times 
with particle's px was discussed. It was noted that due 
to the particular features of the hydrodynamic emis- 
sion function, particles with larger px were, on average, 
emitted earlier. In addition, studies with the "blast- 
wave" parametrization j42| showed that a time differ- 
ence between particles of different masses arises in boost- 
invariant models (the influence of the violation of boost 
invariance on the time and longitudinal shifts has been 
discussed in (T^I). In our case both effects would con- 
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FIG. 5: (Color on-line) Pion-kaon emission functions in 
LCMS for space a) and time b) components, combined into 
the observable asymmetry in PRF c). Red circles are for all 
particles, blue squares for primordial only. 



tribute to the additional asymmetry in PRF: 
Kut = It {rout - Pt^t) , 



(9) 



which would go in the same direction as the spatial one 
coming from radial flow. Even though it cannot be di- 
rectly correlated with the radial flow, it is still very much 
hydrodynamic in nature and we do not consider it as an 
alternative, non-collective explanation of the asymmetry. 

However, the resonance decay process can also intro- 
duce additional time asymmetries, as resonance decays 
occur with a certain time delay and if some particles are 
more abundantly produced by resonances than others, 
asymmetries may arise. In contrast to the effect dis- 
cussed in the previous paragraph, these time delays are 
not "hydrodynamic" in nature and should be treated as 
alternative sources of asymmetry. Moreover, one expects 
pious to be most abundantly produced by resonance de- 
cays, which would mean that they would appear to be, on 
average, produced later than kaons and protons, produc- 
ing asymmetry in the same direction as the flow. Clearly 
the matter requires careful quantitative study. 

In Fig. [5] the time and space emission asymmetries are 
shown on the left and central panels. One sees that both 
are present and are significant. Resonance decays modify 
both distributions and add a complication of long-range 



tails. Therefore, a simple mean of a distribution is no 
longer a good variable to characterize such asymmetries. 
One is forced to use the mean values of the functional 
forms fitted to the distributions around the peak. Tab. lIVI 
shows the summary of these fits. We see that even for the 
primordial particles only, there is already a time differ- 
ence, although smaller than the space one. Introducing 
resonances increases the time asymmetry. Evidently the 
expected effect is seen. Comparing the values in PRF we 
take the difference between "all" and "primordial" cases 
for the time asymmetry as an estimate on how much 
time asymmetry the resonance decays introduce. Both 
for pion-kaon and pion-proton the "non-hydrodynamic" 
asymmetry coming purely from resonance time delay is 
less than 15% of the predicted overall asymmetry. For 
kaon-proton it is less than 25%. Again, the fiow asym- 
metry dominates the calculated asymmetry signal. 



E. Expectations for qualitative trends 

Following the discussion in the previous subsections 
one can formulate several predictions of expected qual- 
itative trends. We will consider pious, kaons and pro- 
tons in specific pT windows, corresponding to the ac- 
ceptance of the STAR experiment at RHIC. This means 
that the velocities of these particles will be fixed, and one 
might assume that the same pious will be correlated with 
kaons for the pion-kaon correlation and with the protons 
for the pion-proton correlation. Since we are correlat- 
ing particles with similar velocity, we will be correlating 
very low-p-r 0.1 GeV/c) pious with medium-p^ kaons 
(« 0.5 GeV/c) and moderate-pT protons (« 1 GeV/c). 
We give a detailed relation between single-particle and 
two-particle sizes in Appendix \^ In terms of the ob- 
servables themselves, that is the two particle variances, 
one expects that Rttk and will be similar and large, 
since they are dominated by the large \ow-pT pion size. 
In contrast Rxp is expected to be significantly smaller. 

The asymmetries show a common feature - with the 
definition of Eqs. (O,® (the lighter particle always taken 
as first) they are all negative - which reflects the fact that 
lighter particles are expected to be emitted closer to the 
center of the system than the heavier ones. 

As for the relations between asymmetries, hydrody- 
namics naturally predicts that the /x^p will be the largest, 
fJ-TrK will be of similar magnitude but smaller, while fj,Kp 
will be much smaller than the other two. It can also be 
shown that the following relation should hold: 



(10) 



The asymmetries for "all" particles for pion-kaon and 
pion-proton pairs should be significantly different than 
for "primordial" only, while for kaon-proton the differ- 
ence should be small. 

In Figs. [HI [71 [S] the predicted system size and emission 
asymmetry for pion-kaon, pion-proton and kaon-proton 
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FIG. 6: (Color on-line) The pion-kaon overall radius a) and 
emission asymmetry b) in PRF as a function of centrality for 
Au+Au collisions at ^/sNM ~ 200 GeV. Close points are for 
all particles, open for primordial only. Circles are same-sign 
pairs, triangles - opposite-sign. 
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FIG. 7: (Color on-line) The pion-proton overall radius a) and 
emission asymmetry b) in PRF as a function of centrality. 
Close points are for all particles, open for primordial only. 
Squares are same-sign pairs, triangles (diamonds) - opposite- 
sign. 



pairs respectively, calculated by Lhyquid-|-THERMINATOR, 
are shown. All the expected trends mentioned in this 
section are confirmed. 



F. Centrality dependence 

We first consider the centrality dependence of the sys- 
tem size. The assumed initial conditions have a clear de- 
pendence: the system size and initial temperature grow 
with Nparu both should result in larger sizes thorough 
of the evolution. The hydrodynamic evolution assumes 
identical equation of state and identical freeze-out tem- 
perature for all centralities; so, there is no reason to 
expect that this dependence will be altered in the final 
state. Also, calculations for identical pions show the same 
trend. Finally, all available experimental data, including 
pion, kaon and proton femtoscopy show the same trend. 
Therefore increase of the system size with Npart is ex- 
pected. 

The asymmetry is the result of a hydrodynamic evolu- 
tion. As already mentioned the parameters of the evolu- 
tion do not change with centrality, only the initial con- 
ditions. Since the chemical properties of matter do not 
depend on centrality either, the resonance decay phase 
is not expected to be very different, giving an enlarge- 
ment of the system size by a constant amount. On the 



other hand, as the collisions become more peripheral, the 
initial overlap region is shrinking, so one expects that 
the long-range resonance corona will become relatively 
more important. Since we know that the overall size of 
the system grows with centrality, the absolute value of 
asymmetry will also grow. The ratio of asymmetry to 
the system size is shown in Fig. |9] and indeed we can see 
that the scaling of /i with holds well for all centrali- 
ties, all pair types, both for "all" and "primordial only" 
particles. 



IV. NON-IDENTICAL PARTICLE 
FEMTOSCOPY FORMALISM 

In femtoscopy one aims to measure the space-time con- 
figuration of the emission process in hadronic collisions 
by analyzing the specific behavior of the two-particle cor- 
relation function. A natural variable versus which this 
correlation is measured is the half of generalized pair 
relative momentum fc*, as opposed to analyses focused 
on event structures, which use e.g. azimuthal angle and 
pseudorapidity differences, k* is calculated in the pair 
rest frame, so it is also the momentum of the "first" par- 
ticle of the pair in the PRF. The particles in the pair are 
different, so it is important to define which one is "first" . 
Later we give the conventions which we use in this paper. 
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FIG. 8: (Color on-line) The kaon-proton overall radius a) and 
emission asymmetry b) in PRF as a function of centrality. 
Close points are for all particles, open for primordial only. 
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FIG. 9: (Color on-line) The emission asymmetry scaled by the 
overall system size in PRF. Closed points are for all particles, 
open for primordial only. Circles are for pion-kaon pairs, 
squares for pion-proton, stars for kaon-proton. 



Femtoscopy also requires a precise knowledge of particle 
type, which means experiments wishing to do such anal- 
ysis must have good particle identification capability. 
The correlation function is defined as: 



C{Pa,Pb) 



P2{Pa,Pb) 
PliPa)Pl{Pb)' 



(11) 



where P2 is a conditional probability to observe a par- 
ticle with momentum pi, if a particle of momentum Pa 
was also observed, while Pi is a simple probability to ob- 
serve a particle with a given momentum. Note that this 



definition is general and applies equally well to the fem- 
toscopic correlation function and to other two particle 
correlations e.g. event structure ones. 

The experimental femtoscopic correlation function is a 
specific form of Eq. (fTTj) . written as: 



C(fc*) = 



J A{pa,pb)d{k* - \{pl-pl))(f'Pad^Pb 
j B{pa,Pb)5{k* - \{pl-pl))<Ppa<Ppb 



A{k* 



B{k*y 



(12) 



where A{pa,Pb) is the distribution of correlated pairs (i.e. 
both particles coming from the same event) of particles of 
type a and h and B{pa,Pb) is the same distribution, but 
particles are not correlated (i.e. come from two different 
events) . Note that the argument of C is changed to half of 
the pair relative momentum k*. Also note that with this 
definition the correlation function will contain not only 
femtoscopic correlations, but also all other event- wide 
correlations projected to the two-particle space. These 
include: elliptic flow V2, global event energy and momen- 
tum conservation, resonance decay correlations (if a and 
b are different and can be products of the decay of a given 
resonance, e.g. products of the A++ resonance in the 7r+ 
- p correlation function) , residual correlations (remnants 
of the femtoscopic correlations between parent particles, 
which decayed weakly into the particles of interest, e.g. 
residual correlations between A and p feeding into the p~p 
correlations), jets, etc. There are numerous experimen- 
talists recipes to construct B in such a way that these 
correlations are included there. In that case, dividing A 
and B also divides out the correlations. One should also 
correct for other non-femtoscopic effects, so that one is 
left with pure femtoscopic correlation in C. The latter 
is desirable, since it is the femtoscopic-only effect that is 
usually computed in models. But such correction pro- 
cedures are never fully effective; so, one must take it 
into account when comparing the "compound" correla- 
tion functions from the experiment with the "pure" fem- 
toscopic CF from models. The details of such proce- 
dures are clearly experiment dependent and are beyond 
the scope of this paper. We only note that some models 
may include the effects of global correlations (e.g. elliptic 
flow, energy- momentum conservation); so, they can be 
used to model such effects. We discuss how this can be 
done later in the chapter. 

In models one defines the correlation function via the 
single and two-particle emission functions: 



Sa{xi,pi) = 

S{xi,Pi,X2,P2, ■.■,XN,PN)dX2dp2—dXNdpN (13) 

Sab{xi,Pi,X2,P2) = 
S{xi,Pi,X2,P2, .■.,XN,PN)dx-idp3...dxNdpN (H) 



which are interpreted as a probability to emit a particle 
(a pair of particles) from a given space-time point with 
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a given momentum. With this definition one might sub- 
stitute Sab for P2 and Sa for Pi in ([TT|) . Such a correla- 
tion function is a 14-dimensional object (7 independent 
components per particle: 4 space-time and 3 momenta). 
Assumptions are needed to reduce the number of dimen- 
sions. We describe them below. 

In principle, the source emission function Sab should 
reflect all the physics aspects of the particle emission 
process, including the proper symmetrization for bosons 
and fermions, as well as the influence of the two-body 
and many-body Final State Interaction (FSI). However, 
commonly used heavy-ion collision models do not include 
these effects. Instead, one assumes that each particle's 
emission process is independent, or in other words that 
an interaction between the two final-state particles after 
they are created is independent from their emission pro- 
cess. The introduction of this factorization of the FSI 
and two-particle wave function symmetrization gives the 
equation: 



C(Pi,P2)= / 5ab(pi,xi,P2,X2) rf'*xi(i''x2(15) 



where 4' is the pair wave function. Particle types A and B 
are known, so the momenta have only three independent 
components, while for positions all four components are 
independent. We also mention that the equation (|15p has 
strong similarities to the Fermi equation used to describe 
the /3 decay process [77| . 



A. Pair wave function 

The pair wave function ^ describes the behavior of a 
pair of particles, one of type A and another of B. In non- 
identical particle correlations we use a particular form of 
^I^, which corresponds to the following physical scenario: 
two particles A and -B, which, shortly after they are pro- 
duced in heavy-ion collision, interact via the FSI, in our 
case Coulomb and/or Strong. After this interaction they 
propagate to the detector as plane-waves. In that case, 
the particular form of is the solution of the scattering 
problem, viewed with the reversed time direction. We 
also use equal time approximation which assumes the 
particles were born at the same time in PRF (see [13, [l3| 
for a detailed description and estimation of the system- 
atic error introduced by such assumption). We factorize 
the wave function into the part describing the motion 
of the pair as a whole (a function of pair total momen- 
tum and "average" emission point), and the component 
describing the interaction itself - dependent on pair rel- 
ative momentum k* and separation r* . The first com- 
ponent produces only an additional phase, which does 
not influence the modulus of the wave function. Since in 
our study we are only interested in the modulus, we can 



neglect this component, and we are left with [K 



-ik*r* 



F{-zr,,l,tC) + fcik*) 



(16) 

where Ac is the Gamow factor, ( — k*r*{l + cos 6**), 
77 — l/{k*ac), F is the confluent hypergeometric func- 
tion, G is the combination of the regular and singular 
s-wave Coulomb functions and fc is the strong scat- 
tering amplitude, modified by the Coulomb interaction. 
0* is the angle between the pair relative momentum k* 
and relative position r* in PRF, while ac is the Bohr 
radius of the pair, equal to 248.52 fm, 222.56 fm and 
83.59 fm for pion-kaon, pion-proton and kaon-proton 
pairs respectively, and is negative for opposite-charged 
pairs. For identical particles 5* must also be properly 
(anti-)symmetrized. This equation is valid in the outer 
regions of the Strong Interaction potential, and neglects 
the components for angular momentum ^ > 1, the latter 
is a valid approximation for small fc*. 

In femtoscopic analysis we assume that we know ^' 
with infinite accuracy, so that we can try to invert 
Eq. to obtain, from the measured correlation func- 
tion C, the information about the emission function Sab- 
In this work we focus on non-identical combinations of 
the most abundant stable hadrons measured in heavy- 
ion collisions: pions, kaons and protons. For each of the 
combinations (pion-kaon, pion-proton and kaon-proton) 
there are four charge combinations: two of the same sign 
and two of the opposite sign. The wave function for 
both same-sign (and both opposite-sign) combinations 
are identical. We note that the wave function charac- 
teristics are indeed well known for all combinations, ex- 
cept for the opposite-sign kaon-proton. The strong in- 
teraction in this system is interesting in its own right, 
and is the focus of intense theoretical investigation (see 
e.g. [Z^)- The femtoscopy technique can be useful in 
this regard - by inverting the problem and assuming that 
we know the source distribution (from other femtoscopic 
measurements) we can invert Eq. (|15p to deduce the pa- 
rameters of ^' from the measured correlation function. 
Similar technique can also be used to study the strong 
interaction in the tt+tt" system precisely. Such studies 
are beyond the scope of this paper (see, e.g., [lOl W^). 

For the systems analyzed in this work, fc can be 
parametrized in the effective range approximation by: 



Ic{k*) 



Jo ^ 



— h{k*ac) 
ac 



ik*Ac{k* 



(17) 

where the /q is 0.137 fm for same-sign pion-kaon pair, 
—0.071 fm for opposite-sign pion-kaon pair, —0.148 fm 
for same-sign pion-proton, 0.112 fm for opposite-sign 
pion-proton, —0.360 fm for same-sign kaon-proton pair. 
The effective radius c?o can be put equal to for all con- 
sidered pairs, at small k* where the I//0 term dominates. 

The full form of the wave function (jl6p must be used 
when calculating correlation functions to be compared 
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with data. However, for theoretical calculations it is 
sometimes instructive (and faster numerically) to neglect 
the strong interaction. This is acceptable since the strong 
interaction is expected to be small for the pairs of inter- 
est, except for opposite-sign kaon-proton pairs, which we 
will not include in the model calculations. For simplicity 
we will use the Coulomb-only wave function in the dis- 
cussion in the next paragraph, but the conclusions hold 
for full Coulomb-|-strong wave function as well. 

For the discussion it is important to provide the form 
of the F function explicitly: 



F{a, 1, z) = 1 + az + a{a + l)z^/2!^ -t- 



(18) 



B. Emission function 



The FSI correlation, described by p6)) . depends only 
on relative momentum k* and separation r* of the pair 
(the angle between the two vectors is 6*). The first sim- 
plification of Sab from is to change to the relative 
variables and integrate out the sum ones: 

SAB{k*,r*) = J J '5'as(xi,Pi,X2,P2) 
x<5(fc*-l/2(pj-p*)) 

xSir*-ixl-x*,)) 

d^xid^X2 (19) 

The integration over space is done over the full variable 
range. On the other hand, particles' momenta are mea- 
surable, so it is possible to define the boundary momenta 
Pmin and Pmax for which the integration is done. In fact, 
for identical particle femtoscopy it is common to define 
several kx = l/2{p\, + p^) ranges and create separate 
correlation functions for each of them. In this way, the 
information about the P dependence of Sab is not com- 
pletely integrated out and can still be inferred. Up to now 
similar binning was not possible for non-identical parti- 
cle correlations at RHIC, because of small statistics and 
limited pt acceptance with good PID coverage. However, 
when one compares the model and the experiment data, 
one must take care to restrict the P integration range 
at least to the one dictated by the pT acceptance of the 
experiment. 

The emission function is a 7-dimensional object. 
We use it to rewrite the equation (IT5|) . putting the specific 
form of the wave function (ITCl) as well: 



C{k*) 



SAB{k*,r* 



{ + ).AB 
-fc* 



(20) 



This form can be explicitly used to calculate the corre- 
lation function from models. However, it is very rare for 
a model to provide the full two-particle emission func- 
tion. The existing models fall into two categories: they 



either provide an analytic form of a single particle emis- 
sion function, or they provide information only about the 
produced particles. 

In the first case we assume that each particle's emis- 
sion process is independent. Then, the two-particle emis- 
sion function Sab can be constructed from single particle 
emission functions via a convolution: 



SABik*,r*) = y S'a(Pi,Xi)S'b(p2,X2) 



xSik* --{pi+ P2))S{r* - xi + xa) 



X d xi d*X2d^pi d^p2 ■ 



(21) 



In case of identical particles {A = B) several simplifi- 
cations can be made. The convolution of two identi- 
cal Gaussians is also a Gaussian with a multiplied by 
\/2. Femtoscopy can provide information about the two- 
particle emission function only, but with the simplifying 
assumption above, the a of the single particle distribution 
can be inferred. For non-identical particles, A B, the 
above simplified method of comparison is not applicable. 
The comparisons with models is more complicated, but 
retrieving single-particle source sizes is still possible, pro- 
vided a complete set of measurements is performed. The 
formula is used in Sec. IVIII and derived in Appendix [X] 
Also, since generally Sa ^ Sb, the Sab can yield a non- 
zero mean value of the separation vector (r*). Later 
in this Section we explain how this mean value can be 
accessed experimentally and argue that this is an im- 
portant and unique piece of information accessible via 
non-identical particle femtoscopy only. 

In the second case, of a model discretely producing 
particles, Eq. (|20l) is evaluated via the Monte-Carlo pro- 
cedure: 



C{k*) 



{+),AB 



(r*) 



(22) 



Note that if the particles from the model are produced in 
a correlated way (e.g. with energy and momentum con- 
servation for the full event, with energy and momentum 
conservation for resonance decay, etc.) these are not de- 
stroyed and are still present in C. In other words, this 
method does not require the simplifying assumption of 
Eq. (I?T|) . In this work we use this method to calculate 
model correlation functions. Its practical implementa- 
tion is described in Sec. IVl 



C. Correlation function and asymmetry signal 



Eq. ([15]) essentially defines the correlation function as 
a pair wave function averaged over the source. Using 
the simplified wave function containing only the Coulomb 
part one can write: 



C{k*) = Aciv)[l + 2 (r*(l -I- cosr)) /a, + 



(23) 
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FIG. 10: (Color on-line) An example of the non-identical par- 
ticle correlation function dependence on the size of the sys- 
tem. Solid lines represent same-sign pion-kaon correlation 
functions for a Gaussian source with sigma from 7.0 fm (low- 
est) to 13 fm (highest) and a shift of —5 fm. The Gamow 
factor is shown as a dashed line for comparison. 



where averaging is done over pairs. In this example we 
illustrate the origins of the asymmetry effect by consid- 
ering same-sign pion-kaon pairs. In such case, Oc is pos- 
itive, Ac is negative and r*(l -I- cos 6**) is by definition 
always positive. For a point source (r* — 0) the overall 
correlation effect R = |C — 1| would be maximum and 
equal to 1 — Ac- For our pairs Ac — 1 is negative, while 
2 (r*(l -I- cos 6'*)) /ttc is positive, so the two compete with 
each other. As the size of the system grows, so does aver- 
age r* , and the correlation effect R decreases. Therefore, 
the correlation function is sensitive to the source size; so, 
we expect to be able to measure the size of the system. 
One has to remember that as the size grows, the CF gets 
less and less sensitive to the system size; therefore, the 
analysis is able to accurately measure only sizes that are 
not too large. To determine whether the "too large" is 
larger than the maximum expected system size at RHIC 
is one of the objectives of this paper. Figure [TU] shows an 
example of how a non-identical particle correlation func- 
tion (in this case same-sign pion-kaon) depends on the 
size of the system. One can see that for a reasonable sys- 
tem size (comparable to the maximum femtoscopic sizes 
obtained in the central Au-(-Au collisions at RHIC) there 
is a noticeable and monotonic dependence on the size, 
indicating that a qualitative femtoscopic analysis should 
be possible. 

As mentioned in the previous paragraph, the two- 
particle emission function can yield a non-zero mean 
value of the separation vector (r*). We now discuss how 
it can be observed in data. From Eq. (|23l) the correla- 
tion function depends on the angle 9*, between relative 
momentum k* and relative position r*. When the two 
are aligned (cos6'* > 0), the correlation effect is smaller, 
than when they are anti-aligned (cos 6* < 0). The former 
configuration means that the particles of the pair, when 
born, immediately start to fly away from each other, so 
their effective interaction time is shorter. In the latter 
case, when they start to fiy towards each other, pass 



close to each other and only later fly away. Angle 9* 
is not accessible experimentally, but it does influence R. 

Particles' momenta, measured experimentally, can be 
used to calculate the relative momentum k* and pair 
total momentum K, corresponding to the velocity v. 
The angle between the two is tp. One divides the ob- 
served pairs in two groups: one having k* and v aligned 
(cosV' > 0) and another, having k* and v anti-aligned 
(cosV' < 0), and creates two correlation functions, (7+ 
and C_ . If in the pair sample used to calculate the C+ we 
have a majority of pairs which also have cos 9* > (and 
in the C_ sample the majority of pairs have cos9* < 0), 
than C+ and C_ would differ. In this particular case C+ 
would show a smaller correlation effect and C_ a larger 
one. When plotting a "double-ratio" C+/C_ one would 
see a signal deviating from unity. For same-sign pairs 
it would go above unity, while for opposite-sign pairs it 
would go below unity. 

If we see a non-zero "double-ratio" , it means that cos 9* 
is somehow correlated with cos 4'- They are connected 
via a third angle of interest: the angle (j) between the 
pair velocity v and the pair relative position r* . When 
we consider only the projections of all these angles on the 
transverse plane, we have trivially: 

i; 9* + (j). (24) 

For the average cosines of these angles, we can write: 

(cosip) = {cos 9* COS (j)) . (25) 

By definition for C+ the left-hand side of Eq. (|^5)) is pos- 
itive. We assumed that C+ shows larger correlation ef- 
fect, so {cos 9*) is negative. The only way that the equa- 
tion can be fulfilled is to also have cos0 < 0. In other 
words it is required that, on average, r* is anti-aligned 
with V. This is a crucial point, so let us restate it. If 
we see a non-unity double ratio we can conclude that 
the average relative position direction is correlated with 
the pair velocity direction. That means that we can ac- 
cess, via a rather straightforward procedure, the mean of 
the two-particle separation distribution, which is allowed 
to be non-zero for non-identical particles. The emission 
asymmetry (r*) is a three- vector, while the considera- 
tion above only mentions a single direction (the direction 
of the pair velocity v). The argument is more general: 
we can replace v by any other direction, defined in the 
Pair Rest Frame, and repeat the argument to obtain the 
same conclusions. The simplest generalization, which we 
have shown to have important physics motivation, is the 
decomposition of v into components: the longitudinal 
"long" (along the beam axis) and transverse "out" (per- 
pendicular to the beam axis). By performing the asym- 
metry analysis versus these two directions we can obtain 
information about the "out" asymmetry fiout = {f*out) 
and the "long" asymmetry /i;ong — (j'long^- ^'^^ com- 
pleteness we also use the third direction: "sidewards" or 
"side" , perpendicular to the other two to study the "side" 
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FIG. 11: An example of the non-identical particle "double- 
ratio" dependence on the emission asymmetry. Solid lines 
represent same-sign pion-kaon "double-ratios" for a Gaussian 
source with sigma of 10.0 fm and the asymmetry of —12.0 fm 
(lowest) to —4 fm (highest). The arrows show the asymptotic 
values (see text for details). 



asymmetry fiside = {^*side)- Each of them carries impor- 
tant physics information or is useful as an experimental 
cross-check; the details will be discussed in the following 
Sections. 

We discussed the possible physical origins of the emis- 
sion point asymmetry and its significance in Section [IIII 
Now we discuss the general properties of the "double- 
ratio" observable. It is allowed to go both above and be- 
low unity, which means either positive or negative emis- 
sion asymmetry (the average emission separation being 
aligned or anti-aligned with the velocity direction). It can 
also be unity, meaning no emission asymmetry. More de- 
tailed analysis also shows that, for a fixed source size, in- 
troducing larger and larger emission asymmetry produces 
larger and larger deviation from unity of the "double- 
ratio" signal. This means that one can determine the 
existence of emission asymmetry and also measure its 
magnitude. The technical details of how it is done are 
presented in the next Section. The illustration of the 
"double-ratio" behavior, for a fixed system size for same- 
sign pion-kaon pair, is shown in Fig. 1111 

It is instructive to derive the equation for asymptotic 
behavior of the double ratio C+ / C_ at k* going to zero ^, 
[§]. We focus on the function F. In this limit, we neglect 
all components with \/ac or \/k* in powers greater than 
1. We have: 



i^=l + 2— + 2f^ + 
ac k*ac 











+ 


_{k*) 







which gives the correlation function: 

C\k' = {AcF*F) ^Ac(l + 2^ + 

V ac k*ac 

One notes that: 



(26) 
(27) 



We now consider C at a fixed k* and cos('0): 
C(/c*,cosV) = ^c f l + 2^ + 2cosV-^^^ ) , (29) 



ac 



ac 



and using the uniformity of cos distribution for uncor- 
related particles at small k* we have: 



' out/ 



C+ ^ Ac C (cos ■ip)d cos ^1 + 2- — - + 

Jo ac ac 

C- ^ Ac [ C (cos ip)d cos tp = 1 + 2 

J -I ac ac 



^1.^0 = 1 + 2^'^' 
G- ac 



out I 



(30) 



+ klrl cos{0* - lb)) = k* cos ^ {r*^t) . 

(28) 



These asymptotic values of the double-ratio are shown 
in Fig. [TT] as arrows. The formula appears to give an 
easy way to extract an emission asymmetry without the 
need for tedious analysis. However, one must remember 
that in the experiment, the lower the fc*, the higher is 
the experimental uncertainty on the data point. This is 
because of statistics, which falls as fc*^, and because ex- 
perimental effects like momentum resolution result in the 
largest systematic uncertainty in these bins. Therefore, 
one has to perform the full analysis of the double ratio 
in a broad range of k* to reliably extract the asymmetry. 

From Eqs. (1501) one concludes that if one restricts the 
integrals in the definitions of C+(C_) to a cos(-0) range 
close to 1(— 1) one will obtain an even larger asymme- 
try signal. However, in the experiment, the price to pay 
is the loss of statistics, and hence the significance of the 
signal. It can be shown that in order to maximize the sig- 
nificance of the signal one should perform the integrals 
over the full range of cos('!/') (thus minimizing the sta- 
tistical error) with the weight cos(^/') [tI, H^. We come 
back to this crucial point in Section FV Al where we discuss 
the spherical harmonics representation of the correlation 
function, which happens to naturally introduce similar 
weighting |83,|8l|. 

We finish this chapter by discussing the conventions 
used in the analysis. It is important to define and con- 
sistently use these conventions in all steps of the analy- 
sis. First, the order of particles in the pair is important 
for the definition of k* and r* because they are defined 
as the momentum and position of the first particle in 
the pair. We adopt a convention that the lighter parti- 
cle in the pair is always taken as first. If both particles 
have equal mass, the positively charged one is taken as 
first. The second convention is the definition of the dou- 
ble ratio, which can either be C+/C_, or C-/C+. We 
chose the former definition. We also note that in the 
spherical harmonics representation there is no need for 
such a convention. With these definitions the following 
general rules hold. Same-charge pairs have correlation 
functions going below unity, opposite-charge pairs above 
unity. That means that for a given source size asymme- 
try, the double ratio for same-sign pairs will be an inverse 
of the opposite-sign double-ratio. Finally, with these def- 
initions, a "double ratio" below (above) unity means that 
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the lighter particle is emitted closer to the center of the 
system and/or later than the heavier one for same-sign 
(opposite-sign) pair. 



function, since we are not interested in non-femtoscopic 
correlations, but rather in physics effects accessible via 
femtoscopy. 



V. PRACTICAL IMPLEMENTATION OF THE 
FORMALISM 

When using a model that produces individual parti- 
cles (such as THERMINATOR), the integration in Eq. ((T5]) 
is performed via the Monte-Carlo method, iterating over 
particle pairs. In that case, the procedure to construct 
the correlation function closely resembles the experimen- 
tal one. It enables the introduction of some experimental 
effects, such as acceptance or momentum resolution, in a 
straightforward way, as opposed to the purely analytical 
models where this is more difficult. Mathematically, the 
procedure is a combination of Eqs. (|12l) and dTSl) : 



C{k*) 



jA{k*,r*) |«'(fc*,r*)|^ 
J B{k*,r*) 



(31) 



Note that the A and B now depend on space-time co- 
ordinates r* as well, because we are using model data 
in which the emission points are known. Since we are 
dealing with pairs of individual particles we employ the 
Monte-Carlo procedure, which replaces the analytic inte- 
gration by summing over pairs: 



Dpairs 



5{k*g - k*) 



(32) 



where the two sums are performed over two sets of pairs 
N and D and the 5 function ensures that only pairs 
with the correct relative momentum are taken. Note 
that there are two possible scenarios, both of which have 
slightly different interpretations and uses. One can per- 
form the calculation using the same pairs for set N and 
D. Eq. ([5^ is then exactly equivalent to dTSl) . but is done 
via Monte-Carlo integration. Another option is to take 
pairs from the same event as sample N and pairs of par- 
ticles from different events as sample D. In that case, the 
correlation function C contains not only the femtoscopic 
effect but also all other event-wide correlations which are 
present in the model, projected to the two-particle space. 
It is therefore very useful for the experimentalists who 
can study the differences between the two to estimate the 
magnitude of the non-femtoscopic effects. The two cases 
have one more important difference. Usually the CF is 
constructed as a ratio of two histograms: N , a signal, 
which is filled with the weight for every pair, and D, 
the background, which is filled with unity for each pair. 
The errorbars on C = N/ D have different meaning in 
the two cases. In the first one it is just the spread of the 
weight in a given sample. In the second case it is a true 
errorbar, comparable to the experimental one, since sam- 
ples TV and D are statistically independent. In this work 
we are using the first way to calculate the correlation 



A. Correlation function representation 

Femtoscopic correlation functions have been repre- 
sented in two main forms: as a 1-dimensional histogram 
with the magnitude of the relative momentum k* or 
q = 2k* on the axis, or as a 3-dimensional histograms 
with fc*„j, and ki^^^ on the axes. For identical pion 
correlations it is also useful to use a 3D histogram with 
the relative momentum components qout, Qside and qiong 
calculated in the LCMS. Note that out, side and long 
decomposition is possible also for pairs of non-identical 
particles; one has to use the generalized four-momentum 
variable q instead of g: q = q — P{qP)/ . The first 
form requires minimal statistics, but only allows to de- 
termine the ID overall source size. The second one allows 
for the determination of sizes in all 3 directions, but re- 
quires significant statistics. Up to now, the non-identical 
correlation function was represented as a set of two 1- 
dimensional histograms - one for C+ and one for C_ with 
respect to the "out" direction. This allowed for the de- 
termination of ID source size and a study of the "double 
ratio" to access asymmetries. However, if one wanted 
to study "double ratios" for other directions ("side" and 
"long"), one needed to create separate sets of correlation 
functions. 

Recently, a more advanced way of representing the cor- 
relation function, the Spherical Harmonics (SH), was pro- 
posed [s^ . It has several important advantages. It en- 
codes the full 3D information of the correlation in a set of 
ID plots. Generally, this does not need to be an advan- 
tage, because a perfect representation of all the features 
of any 3D function requires an infinite set of com- 
ponents (meaning: infinite set of ID histograms). But it 
so happens that the intrinsic symmetries of a pair distri- 
bution in the femtoscopic analysis result in most of the 
l,m components to vanish. It has also been shown that 
out of these that do not vanish, only those with small 
I contain important information, which means that one 
can safely truncate the decomposition at a rather small 
I without the fear of losing any physics information. It 
is as if the spherical harmonics have been specifically de- 
signed to efficiently represent a femtoscopic correlation 
function. 

The spherical harmonics representation has the ad- 
vantages of both the ID correlation function (because 
it requires less statistics) and a 3D one (because it en- 
codes the important part of the 3D information). How- 
ever, first attempts to apply the decomposition methods, 
which were developed for identical particles, to the non- 
identical particle correlations were not successful. Es- 
sentially, one had to first construct the numerator and 
denominator as 3D histograms (usually in |q|, cos 9q and 
(j)) , then divide them and decompose the resulting 3D cor- 
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relation function, negating the low-statistics advantage 
of the ID representation. The decomposition procedure 
also relied heavily on the symmetries present for identical 
correlations. This presented a problem for non-identical 
correlations, because their primary goal is to study the 
emission asymmetry. This breaks one of the symmetries 
present in identical analysis. Also, the single particle 
acceptance of some experiments produced "acceptance 
holes" i.e. regions of empty bins in the 3D function. Also, 
statistics for non-identical pairs were significantly lower 
than for pion-pion correlations, so filling all the bins in a 
3D correlation function with significant number of pairs 
became a challenge. To solve these problems, a new tech- 
nique to represent the correlation function in spherical 
harmonics was developed 'ss"]. Both the numerator and 
denominator are stored directly in spherical harmonics 
(not as 3D histograms), and the procedure to calculate 
the CF directly from them (again, not involving any 3D 
histogram) has been presented [83|. An additional bene- 
fit of the method is that the covariances between all l,m 
components are explicitly taken into account. Unless oth- 
erwise noted, all further correlation functions presented 
in this work are represented using this method. We also 
recommend that experimental groups use this method. 

As was the case for identical particle femtoscopy, the 
spherical harmonics representation turned out to have 
specific synergies with the non-identical femtoscopic cor- 
relation analysis. The important femtoscopic informa- 
tion is contained in only two Z, m components, while the 
other two can be used as additional cross-checks of ex- 
perimental procedures. If one wishes to analyze the full 
3D information, one needs to analyze only two more I, m 
components - the rest can be essentially neglected as they 
should be consistent with zero, or they will not contain 
additional useful information. 

We now investigate the important components of the 
SH decomposition and their sensitivity to femtoscopic 
information. We perform a simple calculation in which 
we assume the source to be a 3D Gaussian in the LCMS, 
having three different sizes R in three directions [out, 
side and long). It also has a non-zero mean value ^out 
in the out direction: 
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l^outY 



side 



ri 



ong 



Hi 



Ide 



7-)2 

^long 



(33) 



With this source function we perform the integration (jlSp 
to calculate the correlation function and present it in 
spherical harmonics representation. We expect that the 
main femtoscopic information is contained in the follow- 
ing components: Cq, 5RCJ, SRCq and 9fiC|. 

To study the sensitivity to the source size we increase 
the overall size, but keep the radii ratios the same. The 
results are plotted in Fig. [T^l All components show some 
sensitivity, but the Cg is affected the most. In the next 
step (shown in Fig. [131) keep the source size the same, 
while we increase the emission asymmetry {flout)- The 
main sensitivity is in the 5RC| component. More im- 
portantly, no asymmetry results in vanishing ^C\, while 
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FIG. 12: (Color on-line) SH components of the correlation 
function for Rout = 8 fm, Rside = Riong = 4 fm, fiout = —4 fm 

(blue dots). Rout = 10 fm, Rside = Rlong = 5 fm, flout = 

— 5/m (red triangles). Rout = 12 fm, Rside = Riong = 6 fm, 
fJ-out ~ —6 fm (green squares). 



increasing the asymmetry increases the signal in 3?C| 
monotonically, approximately linearly. 3?C| is function- 
ally equivalent to the "out" double-ratio. Obtaining a 
quantitative as well as qualitative information about the 
asymmetry should be possible from the analysis of it (in 
correlation with at least Cq where sensitivity is also seen, 
but to a smaller degree). In addition, from Eq. ([5(7|. we 
concluded that in order to maximize the significance of 
the asymmetry signal one should integrate the correlation 
function with the weight equal to cos(V') Remark- 
ably, the definition of 3?C| is essentially: 

^Cliq) = N /c((7,cos(0,),(/),)sin(0,)cos(0,)dr!,, 



(34) 

where N is the normalization factor, while 9q and 4'q 
are the longitudinal and transverse components of the 
4> angle. The SRC^ component happens to be the opti- 
mal way to maximize the transverse components of the 
asymmetry signal, due to its cos{(j>q) weight, while the 
SftCj* maximizes the longitudinal asymmetry signal, due 
to its cos(6'q) weight: 



nC^iq) = N / C{q,COs{eq),^q)cOs{9q)dnq 



(35) 



Once again it appears as if the spherical harmonics were 
specifically designed for the femtoscopic correlation func- 
tion representation. 

Next, in Fig.[T3]we keep the transverse source size and 
the asymmetry the same, while we change the Riong ra- 
dius. Obviously we see a change in the Cq which refiects 
the growth of the overall system size. But the most sen- 
sitive component is iJiCj which carries information about 
the ratio of the transverse to the longitudinal radii, due 
to its cos^ {6q) weighting. 
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FIG. 13: (Color on-line) SH components of the correlation 



function for Rov 



10 fm, Rside = Rlong 



fm. Emis- 



sion asymmetry ^out is changed from fm (blue dots) via 
— 2 fm (red up-triangles) , —4 fm (green squares), —6 fm (yel- 
low down-triangles), up to —8 fm (violet stars). 
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FIG. 14: (Color on-line) SH components of the correlation 

function for Rout = 10 fm, Rside = 4 fm, )J,out = —4 fm. Rlong 

is changed from 4 fm (blue dots) via 6 fm (red up-triangles), 
8 fm (green squares) to 10 fm (yellow down-triangles). 



Finally (see Fig. [T5|) we keep the sum of the transverse 
radii the same, but we change their ratio. The 3fiC| com- 
ponent, with its cos^(0^) weighting, is the most sensitive 
to these changes, while the others remain practically con- 
stant. 

The above calculations show that by analyzing just 
two components of the SH decomposition (C" and 5ftC|) 
one can already perform a meaningful femtoscopic anal- 
ysis and determine the overall source size and emission 
asymmetry. Adding just two more components (3?C2 and 
5RC|), one can also determine source radii in all 3 direc- 
tions, which would normally require a full analysis of a 
3D correlation function. Full 3D analysis does require 
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FIG. 15: (Color on-line) SH components of the correlation 
function for Riong = 8 fm, jiout = —4 fm. Transverse radii 
sum is kept constant: Rout = 7 fm, Rside = 8 fm (blue dots); 
Rout = 8 fm, Rside = 7 fm (red up-triangles); Rout = 9 fm, 
Rside = 6 fm (green squares); Rout = 10 fm, Rside = 5 fm 
(yellow down-triangles); Rout = 11 fm, Rside = 4 fm (violet 
stars). 



larger statistics than the ID size-|- asymmetry one, as the 
sensitivity to the observables in the I = 2 components is 
smaller than in the I = and I — 1 ones. One should 
stress that in the SH representation one is dealing with 
200 data points (50 points per histogram, 4 components) 
versus the 625000 bins in the 3D histogram representa- 
tion. The huge savings in computation time and method 
complexity do not compromise the physics - one is able 
to obtain essentially the same femtoscopic information 
(3 sizes and emission asymmetry). It means that the 3D 
representation is a particularly inefficient way of storing 
the femtoscopic information, while the spherical harmon- 
ics one seems to be perfectly tailored for that task. We 
add that in all theoretical calculations that we have done 
all the other components were either required to vanish 
from symmetry relations, were consistent with zero, or 
contained femtoscopic information that was already ac- 
cessible via the four main components. In experiment one 
should, in addition, look at the SRCj* (equivalent to the 
"long" double ratio) and 3Cj (equivalent to the "side" 
double-ratio) . The former is expected to show zero asym- 
metry, the latter is required to vanish due to symmetry 
reasons, but their deviations from zero may signal exper- 
imental reconstruction problems. 



B. Extracting qualitative information 

In the previous paragraph the sensitivity of the non- 
identical particle correlation function to the source 
size and asymmetry was illustrated. Dependencies in 
Figs. [121 [131 [m HSl show that source size parameters 
influence all SH components at the same time, and it is 
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a, [fm] 

FIG. 16: (Color on-line) An example of the 2D ^ map ob- 
tained by fitting a non-identical particle correlation function. 
Underlying colored histogram represents the "mesh" obtained 
in the fitting procedure. The lines represent a fitted 2D 
parabola. The point represents the minimum of the parabola 
- the best fit value. 

not possible to analyze them separately. In identical par- 
ticle femtoscopy obtaining Gaussian source size param- 
eters is straightforward: the integral in Eq. (|T5]) can be 
performed analytically if one assumes that S'(r) is a 3D 
Gaussian and does not depend on pair momentum (al- 
though the pairs transverse momentum Kt dependence 
can be recovered by Kt binning; the rapidity binning is 
also possible). To take into account the Coulomb inter- 
action, an approximate wave-function is used in which 
the Coulomb part is factorized out and replaced by an 
averaged function that only depends on the magnitude 
of k* . This procedure [8J|, enables to write a simple 
analytic formula, which can be directly fitted to the 3D 
correlation function and provides femtoscopic radii. 

In non-identical particle femtoscopy the Coulomb in- 
teraction is the source of the correlation, so it cannot be 
factorized out. Performing the integral in Eq. (1151) ana- 
lytically becomes impossible, especially when one needs 
to consider the strong interaction as well. The proce- 
dure must be carried out numerically. One starts with 
the same assumption as for identical particle femtoscopy. 
The source is a 3D Gaussian ([33l) . where the additional 
modification allowing for a non-zero shift in the "out" 
direction is introduced. One then assumes a certain set 
of source parameters (i?o«t, Rside, Riong , fJ-out) and calcu- 
lates the corresponding correlation function, according to 
Eq. (lisp and using the corresponding wave function (|16p . 
As the source function is assumed to be momentum inde- 
pendent, one needs an input momentum distribution (the 
assumption of momentum independence of the source 
function is not required provided the selected {pT,y) is 
sufficiently narrow). This can be achieved by e.g. taking 
particles' momenta from real pairs from data and assign- 
ing the separation r* randomly generated from ([551) . The 
resulting "model" CF is compared to the "measured" 
one via a test. The procedure is repeated for sev- 
eral sets of source parameters to find the one that fits 
the "data" best. This set is taken as the result of the 
fit. The procedure can be refined by ensuring that the 



"model" calculation is done for points which form a reg- 
ular "mesh" in parameters' space. In the simplest one 
parameter fit, one obtains the fit value and error from 
the location of minimum on the parabola. When 2 in- 
dependent parameters are fitted (the most common case 
for non-identical correlations, when we usually fit overall 
radius R and emission asymmetry fi) one creates a 2D 
"mesh" , to which one can fit a 2D parabola - obtaining 
the best fit parameters, their errorbars as well as the co- 
variance between them. An example of such analysis is 
shown in Fig. 1161 One can also employ a minimization 
package (e.g. Minuit) to perform the fitting process. The 
number of points in the CF "mesh" can be kept at mini- 
mum by using a simple linear or quadratic interpolation 
between points [l^, [s^ . 

For the results presented in this work we used a soft- 
ware package CorrFit to perform the numerical fitting 
procedure described above. It was developed for the anal- 
ysis of non-identical particle correlations in the STAR ex- 
periment [3^ . It allows for significant flexibility when 
choosing various parameters of the fitting process. For 
this work we have made the following choices: (1) The in- 
put data were correlation functions in spherical harmon- 
ics representation. Only Cq and SftC^ components were 
fitted; the covariance between the two was taken into ac- 
count. (2) The particular form of the S from Eq. ([55)1 
was taken with the additional constraints that Rside was 
assumed to be equal to Rout, while Riong = l-3i?otit, fol- 
lowing the identic al p article 3D femtoscopic results for 
pions from RHIC [861 ]. Therefore, only two independent 
parameters were fitted: Rout and ^lout- (3) The input 
momentum distributions were taken from the STAR ex- 
periment, ensuring that the momentum acceptance was 
the same. (4) The standard package from R. Lednicky 
was used to calculate the pair weig hts [l3,[il. The input 
theoretical correlation functions were calculated taking 
into account Coulomb interaction only. In such case only 
the Coulomb part of the pair weight was calculated in the 
fitting procedure. (5) One had the possibility to intro- 
duce "pair purity" and momentum resolution corrections. 
Details of the purity correction will be discussed later in 
this Section. The momentum resolution correction was 
not necessary for model studies shown in this paper. 



C. Influence of pair purity 

For femtoscopic analysis one of the main experimen- 
tal issues is the "purity" P of the analyzed pair sample, 
that is the fraction of pairs in the sample that should 
be treated as "femtoscopically correlated". There are 
several reasons why a pair of particles should not be 
treated as correlated. From the experimental side, it 
may happen that one (or both) of the particles in the 
pair has been misidentified (this is why P is traditionally 
called "purity" ) . Another common scenario is when one 
of the particles is a product of a weak decay. In that 
case, it is the "parent" particle that should be treated 
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as "femtoscopically correlated" , leading in some cases to 
rather complicated cases of "residual correlations" , that 
is feed-down of femtoscopic correlation between parent 
particles into the daughters' correlation function. Re- 
cently, such cases have been studied in detail in baryon- 
baryon correlations [lllls^l- However, in systems consid- 
ered in this paper, containing a meson, residual correla- 
tions are smeared away by the decay momenta, so such 
pairs should be treated as "not femtoscopically corre- 
lated" . These experimental components of P are clearly 
experiment-dependent and it is up to the experiment to 
correct for such effects (or at least provide a realistic es- 
timate of P) . We will not address it further in this work. 

However, another contribution to P remains, the esti- 
mate of which is model dependent. If one of the particles 
in the pair conies from a strongly decaying resonances 
that lives very long, e.g. the w meson, we may need to 
treat it as "not femtoscopically correlated". That is be- 
cause the source of the correlation, the pair wave function 
^' is usually peaked at low r* values, while at large val- 
ues it produces no correlation. This can be dealt with 
in two ways. First is to assume, in the data analysis, a 
source function that perfectly describes such long-range 
r* tails. However, this is difficult to do, and it is very 
dependent on the model that we use to model such tails. 
Another solution is to treat the particles in the tails as 
"not femtoscopically correlated" - in other words to pro- 
vide a model estimate of the decrease of P coming from 
long-lived resonances. This is usually done by assuming 
that a source is a 3D sphere with a Gaussian profile, and 
counting any pair which is outside this sphere as non- 
correlated. We note that the "non-gaussian effects" have 
been seen in all RHIC experiments in the identical par- 
ticle correlations [s^, [s^ . They have been attributed, at 
least in part, to long-lived resonances [111, and recently, 
new techniques have been proposed to analyze them in 
detail [89*1. Such analysis are not yet possible in non- 
identical particle analysis, so we limit ourselves to the 
simple model estimation of P. 

We characterize the overall "purity" of a pair sample 
by a percentage of pure pairs P. It can be (and in the 
experiment it usually is) a function of fe*. Assuming that 
the "non-pure" pairs are not correlated, the correlation 
function can be trivially corrected for purity: 

Ccorrected(k ) = {Cmeasured(k ) 1)/P{k ) + 1. (36) 

Employing the formula (|36l) requires a precise knowledge 
of the fraction of correlated pairs P in the measured sam- 
ple. Any uncertainty in its value will be a source of the 
systematic error. In contrast to the identical particle 
femtoscopy, P cannot be easily inferred from the fit, in- 
dependently of the source parameters a and /i. By in- 
specting Fig. [TU] and Eq. ([55)1 one concludes that lowering 
purity influences the correlation function in a way that is 
similar to the changes introduced by varying the source 
size. One can obtain a satisfactory fit to the correlation 
function by adjusting either one of these parameters, so 
treating them both as free makes them highly correlated. 
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FIG. 17: Values of the correlation function fit (size a) and 
shift b) ), as a function of a pre-defined fraction of primary 
pairs. 

While possible, the independent determination of a and 
P requires significant statistics and very good control of 
the non- femtoscopic background p^ . [35| . 

Usually, the practical solution is to come up with the 
best estimate of P by other means (both experimental 
and theoretical), and fix its value for the fitting. By 
performing the procedure for several fixed values of P, 
within the reasonable uncertainty range, one obtains the 
systematic error on the fit values coming from the purity 
estimate. As an example on Fig. [T7] we show how the fit 
values for the same input correlation function change as 
one adjusts the fixed pair purity P. Within the reason- 
able range of P variations of the order of 20% the depen- 
dence of a and /i on P is monotonic and noticeable. We 
expect that the purity estimate will be the major source 
of the systematic error. We will address the theoretical 
part of the estimates of pair purity in Section IVI Bl 

D. Influence of the momentum resolution 

Momentum resolution will influence the femtoscopic 
correlation function in a well defined way. Since recon- 
structed momenta will differ from the true ones, the cor- 
relation effect (visible as either a "peak" or a "depres- 
sion" at low k*) will be reduced. The Coulomb correla- 
tion is rather sharp at low relative momentum, so one ex- 
pects the small k* points to be influenced the most. The 
"double-ratio" will be affected as well. One can study 
the influence of the momentum resolution by performing 
a theoretical calculation in the following way. One cal- 
culates C according to ([5^ . but calculates two sets of 
relative momenta k* for each pair. The flrst one /cj^^g 
is used to calculate the second one A:*,^^^^^^ is used 
to determine the correlation function bin in which the 
pair is stored. The second set is calculated by smearing 
the particles' momenta with the parametrization of the 
momentum resolution obtained from the experiment. By 
comparing the "unsmeared" and "smeared" correlation 
functions one can judge the importance of the momen- 
tum resolution effect. One can also employ the same 
technique when calculating the "theoretical" correlation 
functions during the fltting procedure described in Sec- 
tion lVB] In this way the fitting procedure automatically 
corrects for momentum resolution. The precise determi- 
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nation of the momentum resolution parameters is very 
dependent on the experimental features. We will not 
discuss it further in this work (see, e.g. [s^). 



VI. SYSTEMATIC CHECKS OF THE METHOD 

In the previous sections we have made qualitative 
claims that the analysis of the non-identical particle cor- 
relations should enable the estimation of the source size 
and emission asymmetry for pions, kaons and protons 
emitted in the heavy-ion collision. In this section we aim 
to show the quantitative checks of the method. We fo- 
cus on answering two questions: Can the method reliably 
and quantitatively recover the source size and emission 
asymmetry? If yes, what systematic uncertainties in this 
estimation come from the method itself? 

We emphasize that in this Section we focus only on the 
technical aspect of the non-identical particle correlations 
method. So, the THERMINATOR calculations are used only 
as "test samples". 



A. Details of the procedure 




FIG. 18: Two-particle emission function for pion-kaon pairs 
at centrality 10-20%. r*„t in a), r^^^ in b), r;*„„g in c), ri„„ 
in d). Lines are fits to the distributions, see text for details. 



than the 3D identical pion analysis. Instead, we chose to 
focus on asymmetry, which is a unique observable acces- 
sible only via non-identical particle femtoscopy. 



We chose the following procedure: we simulate heavy- 
ion Au-I-Au collisions with the THERMINATOR model. We 
do it for 6 centrality bins: 0-5%, 5-10%, 10-20%, 20- 
30%, 30-40% and 40-50%. We calculate the full set of 
non-identical correlation functions: same and opposite 
sign pion-kaon, same and opposite-sign pion-proton and 
same sign kaon-proton for these centralities, producing 
6x5=30 independent correlation functions. Then, we pro- 
ceed to treat these functions as if they were experimen- 
tal ones: applying the purity correction and fitting them 
with the CorrFit software. In the end we compare the 
input values from the model with the fit results and see 
if they match. 

We simulated 50K events for each centrality. The pa- 
rameters of the model were the same as the ones used 
in [26l |50| . which are known to reproduce both single 
particle spectra and overall femtoscopic sizes from iden- 
tical particle correlations. The correlation functions were 
calculated by combining particles from these events into 
pairs and employing Eq. (I32p . Spherical harmonics rep- 
resentation of the CF was used. Then, each function 
was fitted with the CorrFit program, assuming that the 
two-particle source is a 3D Gaussian in LCMS, accord- 
ing to Eq. The sideward and longitudinal sizes of 
the source were fixed to be equal to Rout and l.iRout 
respectively, according to the results obtained for identi- 
cal pions. In this way the fitting procedure had only two 
parameters: overall size Rout and emission asymmetry 
A'otit- We note that while the full 3D analysis, with Rside 
and Riong as independent free parameters of the fit, is in 
principle possible, we do not discuss it in this paper, as 
we consider it less interesting: it would provide 3D infor- 
mation about the source, but with much lower precision 



B. Characterizing model input 

Let us inspect a typical two-particle emission func- 
tion Sab coming from the THERMINATOR model, shown 
in Fig. 1181 The distributions in out, side and long direc- 
tions show different behavior. Only the side one is well 
described by a Gaussian; attempts to fit the out and long 
by Gaussians are shown as dashed lines. The long can 
be described by an exponential hyperbola function: 



Jeh = exp 



long 
long 



(37) 



where a and a are parameters. For the description of the 
out one needs an asymmetric exponential hyperbola: 



fas 
I EH 



exp 




fJ-outY 



+ {^ + Cir*out-f^out)), 

(38) 

where aout, fJ'out, ciout and C are parameters. We have 
confirmed that these functional forms are general, that 
is, one is able to fit them to the emission functions at 
all centralities and for all pair types. In contrast, the 
fitting procedure that we have just described in the pre- 
vious Section assumes a proper Gaussian source distri- 
bution. Therefore, one needs to find a suitable variable 
which will enable the comparison between the model in- 
put, which is non-Gaussian, with the fit output, which is 
postulated to be a Gaussian. In principle, one can think 
of removing the Gaussian assumption and actually us- 
ing the functional forms mentioned above for the fitting. 
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This is certainly possible, both from the point of view 
of the technique as well as the software tools. However, 
it brings two potential problems. First of all, while the 
functional forms are able to fit the THERMINATOR data 
very well, there is no reason to assume it will be the 
case for other models. So, by choosing these particular 
forms one would introduce a strong model dependence. 
Moreover, it is not known how the parameters of these 
functions could be compared to the Gaussian radii from 
identical particle femtoscopy. In this work we access the 
overall directionally-averaged size of the system only, so 
the details of the source function dependence in the three 
directions separately will be lost. In this work we decided 
to use the Gaussian assumption, having in mind that it 
may introduce systematic effects, when comparing the fit 
values to input two-particle distributions. 

We are then presented with the following situation: 
the "experimental" fit procedure produces the direction- 
averaged source size R and the emission asymmetry in 
out direction fiout- On the other hand, to fully describe 
the theoretical model input, one needs the following pa- 
rameters: aout, aout, Ci Mout, CTside, Tlong and along- To 

compare the two, one needs to find common variables 
which can be compared, and which will tell whether our 
"experimental fit" reproduces the "model input" . 

In traditional HBT one compares 1-dimensional fenito- 
scopic radii, defined as a ti of a ID-Gaussian approxima- 
tion of the single-particle emission function (usually in 
PRF). We wish to define an analogous variable for non- 
identical study, but defined in LCMS to facilitate the 
comparisons to the 3D identical particle interferometry 
results. The task is complicated by the fact, that neither 
the "model input" nor the "experimental fit" produce a 
source which is a perfect Gaussian. Nevertheless, one can 
plot the source distribution dN/dr in both cases. Then, 
one fits the distribution with the Gaussian formula, mul- 
tiplied by the proper Jacobian: 

FG(r)«cxp(--^)r2. (39) 

As will be discussed later, one must restrict the range 
of this fit to the low-r values, as these are the ones that 
contribute to the femtoscopic effect. The contributions 
from the large-r part of the source must be dealt with 
separately. The lower-right part of Fig. [TH] shows the 
example of the model ri„„ distribution; the line is a fit 
according to (l39l) . done in the range — 20 fm. Since both 
Tout and riong distributions are manifestly non-Gaussian, 
so is the Tmv distribution"'. Therefore, the fitted i?™^"* 
is only characterizing the general size of the system. The 
same fitting procedure can be applied to the rinv dis- 
tribution produced by the "experimental fit" , producing 



* It can be shown that ri„„ distribution is a Gaussian only if all 
three rout, Tgids and riong are Gaussians and additionally all 
three have the same width. 
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FIG. 19: (Color on-line) The correlation between input val- 
ues obtained from model and fit values. Panel a) shows 
overall system size, b) - emission asymmetry in the out di- 
rection. Open points show calculations for primary parti- 
cles only, closed points - for all particles. Violet circles are 
same-sign pion-kaon, red up-triangle: opposite-sign pion- 
kaon, green squares: same-sign pion-proton, yellow down- 
triangles: opposite-sign pion-proton, blue stars: same-sign 
kaon-proton. 




FIG. 20: An illustration of the long-range tails produced by 
THERMINATOR. Points: separation distributions for pion-kaon 
pairs, a): rout, b): Tside, c): riong, d) ri„„. Lines are fits to 
the region near the peak, using the functional forms: f^^'u for 
out, Gaussian for side, /eh for long and fc for inv. 



the i?^** "experimental" value. For a more detailed dis- 
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TABLE V: Estimate of the theoretical purity P for ah pair 
types, versus centrality. See text for details. 



E 



ttK 

Tip 

Kp 



0-5 5-10 10-20 20-30 30-40 40-50 



0.83 0.82 
0.80 0.79 
0.96 0.95 



0.80 
0.77 
0.95 



0.76 
0.75 
0.94 



0.73 
0.73 
0.93 



0.70 
0.71 
0.93 



cussion of the approximate relations between ID and 3D 
source size parameters, the relations between values in 
PRF and LCMS as well as a more detailed description of 
estimating Rav we refer the reader to Appendix [X] 

One also needs to define, what does "emission asymme- 
try" ^out mean in the case of the model. One can come 
up with several definitions: (a) the "mean" of the r*„j 
distribution, (b) the "mean" of the r*^^^ distribution, but 
only taking into account parts of the distribution near to 
the peak (neglecting the long-range tail), (c) the position 
of the maximum in the r*y^^ distribution. All three are 
correlated, but (a) shows large sensitivity to long-range 
tails and therefore is not well defined. Out of the other 
two (b) should be exactly the same as (r*^^), which is 
the value used in the theoretical formalism, while the re- 
lation of (c) to is less straightforward. So, we chose 
to use the definition (b). 

In Fig. [inithe comparison between "input" and "fitted" 
radii is shown, for all considered pairs, both for primor- 
dial and all particles. A clear correlation is seen between 
the two. However, most of the results do not lie on the 
"perfect" x — y line. Within each pair system the corre- 
lation is clear and monotonic. Moreover, the deviations 
from the "perfect" curve seem to be a systematic shift, 
similar for all points of the same pair type. The shift 
from perfect values for a given pair type is seen both for 
the Rinv and the iiout variables. 

One obvious candidate for the explanation of this shift 
is pair purity P. The fit procedure assumes that all the 
pairs in the source come from the region with Gaussian 
density profile. On the other hand, a realistic model, such 
as THERMINATDR, clearly shows that there are significant 
non-Gaussian long-range tails in the separation distribu- 
tions; this is illustrated in Fig. [20l One can perform a 
simple calculation. Take the model radius R^^ obtained 
from fitting the model separation distribution with p9p 
and by integrating this distribution one obtains the num- 
ber of pairs Nq within the Gaussian core of the source. 
No is then the area below the fit curve on the lower-right 
panel of Fig. [501 This can be compared to the total num- 
ber of pairs Na obtained by simply counting the number 
of model pairs, or in other words - the area below the 
points on the lower-right panel of Fig. [501 The value 
fpc — Ng/Na is the number of "femtoscopically corre- 
lated" pairs and should be simply treated as "purity" in 
the sense discussed in Section IV CI One needs to cor- 
rect for these "non-gaussian" effects. In order to do that 
one calculates fpc and then treats it as a fixed parame- 
ter in the fitting process. Note that the fpc estimation 
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FIG. 21: (Color on-line) Same as Fig. ll9l but the fit procedure 
is done with adjusted purity values (see text for details). 



procedure by design takes into account two effects, which 
cannot be easily disentangled - the pairs in the long-range 
tails and the fact that the shape of the core system is not 
a Gaussian in a ID representation. We have performed 
the estimation, based on the input model pair separa- 
tion distributions and find that fpc depends both on the 
pair type and centrality. The exact values are given in 
Tab.|Vl These values are obviously model dependent and 
are strongly influenced by the size of the long-range tails 
in the separation distributions. It is well known that the 
source of such tails are, to a large degree, particles com- 
ing from the strongly decaying resonances. It is therefore 
important to use a model like THERMINATDR, which fully 
includes all known resonance propagation and decay, to 
determine purity. A clear trend exists for all pair types: 
the fpc get smaller as collisions get more peripheral. 
The long-range resonance tails and non-Gaussian effects 
get relatively more important as one moves away from 
central collisions. 

The "experimental fits" have been redone, this time 
with "purity" fixed to the listed fpc values. The re- 
sults are shown in Fig. [51] All results now lie close 
to the "perfect" x = y line within the statistical er- 
ror. This shows that, having in mind the caveat men- 
tioned above, the analysis technique is able to produce 
reliable results, and that the results of the "experimental 
fitting" procedure do indeed provide a valuable informa- 
tion about the particle-emitting source. In particular, 
both the direction-averaged size of the source and the 
emission asymmetry can be reliably recovered. 
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FIG. 22: (Color on-line) Fitted parameters for all pair types, 
versus centrality. System size (a) in panel a), emission asym- 
metry (fi) in b). Colors and symbols are the same as in Fig. ll9l 
The bands represent the change in the fit parameters as the 
fixed purity P is changed by ±0.05 for pion-kaon and pion- 
proton pairs and ±0.03 for kaon-proton. 



The inherent uncertainty in the determination of fpc 
is the source of the systematic error introduced by the 
method itself. In Fig. [52] we show how the fit value 
changes when purity is varied within a reasonable range 
(±0.05 for pion-kaon and pion-proton and ±0.03 for 
kaon-proton). This results in the following systematic 
uncertainties on the obtained values: for same-sign pion- 
kaon 7% on a and 7% on /z, for opposite-sign pion-kaon 
7% on a and 4% on fi, for same-sign pion-proton 8% on 
a and 6% on /i, for opposite-sign pion-proton 7% on a 
and 4% on fi, for same-sign kaon-proton 5% on a and 
4% on fj.. 

We note that the fpc estimation presented here is 
specific to THERMINATOR and represents the uncertainty 
within this model itself. In addition, the absolute value 
as well as an estimated uncertainty of fpc can be differ- 
ent in other models, e.g. in rescattering codes. The total 
systematic uncertainty of the experimental measurement 
should take this model dependence into account. In ad- 
dition, the experimental "purity" , that is the efficiency 
of particle identification, will also contribute to the same 
uncertainty. 



VII. RESULTS OF "EXPERIMENTAL-LIKE" 
ANALYSIS 

In Section |V] we have described, in detail, the exper- 
imental procedure to analyze the non-identical particle 
correlations. For the model analysis of the THERMINATOR 
output, we have used the complete two-particle method 
for calculating the correlation function and we have ob- 
tained the pion-kaon, pion-proton and kaon-proton cor- 
relation functions. They are, from a point of view of 
a formalism and our analysis methods, identical to the 
correlation functions that one might obtain in the ex- 
periment. We calculated sets of ID histograms which 
correspond to theses functions represented in spherical 
harmonics. We neglected the strong interaction compo- 
nent of the pair wave-function for simplicity and speed 
of calculation. Then we switched this effect off in the 
fitting procedure as well. Obviously, when fitting the 
true experimental functions, one will use the full pair 
wave-function calculation. The strong interaction effect 
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FIG. 23: (Color on-line) Example of the correlation func- 
tions calculated from THERMINATOR±Lhyquid model for cen- 
tral Au±Au collisions. Panels a), c), e) show the Cq compo- 
nents, panels b), d), f) show KCi. Closed circles are same- 
sign pion-kaon, open circles opposite-sign pion-kaon, closed 
squares are same-sign pion-proton, open squares are opposite- 
sign pion-proton, triangles are same-sign kaon-proton. Lines 
show the best-fit "model" correlation functions. 



is known to be small compared to the Coulomb for the 
pairs of interest, so we do not expect any systematic ef- 
fect on the fit values coming from this simplification. In 
addition, in the real experimental correlation function, 
one expects some non-femtoscopic effects in addition to 
the pure femtoscopic one. These need to be dealt with 
on a case- by-case basis, as they will strongly depend on 
experimental conditions. For this discussion we assume 
that any such effects can be identified and that the exper- 
imental correlation function can be properly corrected, so 
that only the femtoscopic effect remains. 

We proceeded to treat the calculated correlation func- 
tions as if they were coming from the experiment. No 
other information, except of the pair purities is used 
in the procedure. For purity correction we used val- 
ues listed in Tab. El as we would have done in the real 
data analysis. The goal of the exercise it to confirm 
that the obtained "experimental fit values" actually cor- 
respond to the true values, shown in the previous Sec- 
tion and obtained directly from the emission functions. 
The example correlation functions together with the fit- 
ted "model" ones are shown in Fig. [531 For pion-kaon 
and pion-proton one can see the positive correlation ef- 
fect for opposite-sign pairs and a negative one for same- 
sign pairs. The 3?C| components also show mirror effect. 
For kaon-proton the correlation effects starts at larger 
k* (is wider) and is more pronounced, as expected from 
a smaller Bohr radius for this pair. The lines show the 
functions fitted with the CorrFit program. 

The fitting described in Sec. IIVI assumes that the 
source is a 3D Gaussian in LCMS. The direct output 
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FIG. 24: (Color on-line) The source parameters in PRF for 
pion-kaon pairs obtained from the "experimental" fitting pro- 
cedure. Panel a) shows source size, b): emission asymmetry. 
Circles show same-sign pairs, triangles: opposite sign. Open 
star is a STAR measurement [s^ at ^/snn = 130 GeV, the 
line represents statistical-|-systematic error (see text for de- 
tails). 



of the fitting procedure is a size of the system in the 
out direction a-^ and the emission asymmetry fi-^ in that 
direction. From these values one can calculate, via the 
relations specified in Appendix \^ all the other source 
size characteristics. In the fitting procedure each pair 
is treated individually, so its velocity is known. One 
can therefore directly determine the direction averaged 
source size in PRF - there is no need to use averaged 
pair velocity in transformation from LCMS to PRF, like 
in Eq. (PTs]) . In Figs. [21 [2S1 [Ml we plot this directionally 
averaged radii in PRF i?*^. Out of the three equivalent 
values: cr/, R^^ and R*^ one must choose one for presen- 
tation of the results. When we present the results of our 
calculations, we chose as it is more natural than cr/ 
(which requires the knowledge of other directions multi- 
pliers in order to be meaningful) and describes the emis- 
sion function in PRF, so there is no need for assumptions 
required in the LCMS to PRF transformation. On the 
other hand, the 3D identical particle femtoscopy analysis 
produces source sizes directly in LCMS, so if we want to 
compare to these results, we will use 

In Fig. [24] we show the "experimental-like" fit re- 
sults for pion-kaon pairs. This should be compared 
to the "true" values from Fig. [Bl The correlation be- 
tween "true" and "fitted" values is also shown directly in 
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FIG. 25: (Color on-line) The source parameters in PRF for 
pion-proton pairs obtained from the "experimental" fitting 
procedure. Panel a): source size, b): emission asymmetry. 
Squares show same-sign pairs, triangles: opposite sign. 



Fig. [211 We can see that even with all the assumptions 
and simplifications that are used in the fitting procedure 
one is able to recover the true system size. The accuracy 
for same-sign pion-kaon function is 8%, for opposite-sign 
is 5%, with the largest deviation of 10%. These are com- 
parable to the statistical error of the fit. For the asym- 
metry the "input" and "fit" values were identical within 
10%. All results were in agreement within statistical er- 
ror of the fit. 

A single measurement is available for the pion- 
kaon system at RHIC, done by STAR at ^/SNN = 
130 GeV nil. To compare it with the values presented 
in Fig. [24l we needed to account for two effects. The 
measurement was corrected for a fraction of non-primary 
particles (particles not coming from the primary ver- 
tex), which is an experimental correction. But, it was 
not corrected for non- correlated primary pairs, a correc- 
tion which we described in this work. Taking the de- 
pendency in Fig. [T71 as a guideline and the purity esti- 
mate of 0.85 for pion-kaon in central Au-|-Au we have 
scaled the reported system size and asymmetry accord- 
ingly. Moreover, the fitting performed in 39] assumed a 
specific shape (size equal in out, side and long directions 
in PRF), from which we recalculated Ravi to be com- 
pared with results in this work. The measurement was 
done at a colliding energy different than the one consid- 
ered in this work, so we plot the result at a correspond- 
ing Npart- One can see that the system size is in perfect 
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FIG. 26: (Color on-line) The source parameters in PRF for 
kaon-proton pairs obtained from the "experimental" fitting 
procedure. Panel a) shows source size, b): emission asymme- 
try. 



agreement with the model predictions in this work; the 
asymmetry seems to be slightly smaller than the pre- 
diction, but in agreement within statistical-|-systematic 
error. 

In Fig. [25] the same results are shown for pion-proton 
pairs. The accuracy of the system size determination is 
5%, comparable to the statistical fit error. The asymme- 
try is reproduced with worse quality, with discrepancies 
of up to 12%, comparable to the statistical error. 

In Fig.[26]the results are shown for kaon-proton pairs. 
The system size determination is good: 3%, better than 
the statistical fit error. The asymmetry is reproduced 
with discrepancies of up to 15%; however, since the ab- 
solute values are small (compared to the pion-kaon and 
pion-proton case), the absolute value of the difference 
is comparable with the systematic error of the fit. We 
conclude that the system size and the emission asymme- 
try can be reliably recovered, with the systematic error 
due to the procedure itself less than 10%. However, one 
must correctly determine the systematic error due to the 
pair purity estimation as well. We note that the tests of 
the method have been performed for Therminator model 
only, so they are not necessarily general. However the 
tests were done for many centralities (system sizes) and 
pair types and the method was found to be working in 
all cases, so we have some confidence that it should work 
for other models as well. One possibility for further stud- 
ies would be to perform similar tests with a model that 
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FIG. 27: (Color on-line) Comparison of overall size obtained 
from the non-identical particle calculation with the RHIC 
data for central Au-|-Au collisions at 200 GeV/c. Closed cir- 
cles are (from left to right) single particle pion, kaon and pro- 
ton radii inferred from non-id simulations (centrality 5-10%); 
open triangles are pion radii from STAR [93| (centrality 0- 
5%), open squares are pion radii from PHENIX [9ll] (central- 
ity 0-30%), open diamonds are proton results from STAR [stI ] 
(centrality 0-10%). 



introduces particle rescatterings. 

One would also like to compare the results of non- 
identical particle correlations with the wealth of data 
coming from femtoscopic analysis at RHIC. This presents 
a complication, since the identical-particle femtoscopy is 
usually presented in 3D form as single particle "HBT 
radii" Rout, Rside and Riong, while for analysis in this 
work we only have averaged two-particle source size 
and emission asymmetry ii^ut- From its definition we 
conclude that e.g. tr^^ should correspond to the vari- 
ance of convolution of two Gaussians: one with variance 
of Rout for pions and the other with variance of Rout 
for kaons at the same velocity. Both the identical par- 
ticle i?'s and the non-identical are defined in LCMS. 
Again we refer the reader to Appendix [X] for explicit re- 
lations between the two. From the input data (either 
Rout, Rside and Riong in LCMS for identical particles 
or (7/ and /i/ for non-identical) we calculate the direc- 
tionally averaged source size Rav By comparing these 
values, defined in LCMS, as opposed to the defined 
in PRF, we avoid the unnecessary approximation com- 
ing from the determination of the averaged pair velocity, 
needed for the determination of (7). Such comparison is 
made in Fig. [27l where data from most central Au-|-Au 
collisions are shown. The open points are experimental 
data for pions, kaons and protons. The closed points 
are single-particle source sizes inferred from fits to the 
non-identical particle correlation functions calculated for 
THERMINATOR. One can see a very good agreement be- 
tween the model prediction and the data. All datapoints 
follow the "my scaling" trend predicted by hydrodynam- 
ics. A direct comparison of non-idnentical particle sizes 
and asymmetries between model and data would be an 
even better test. However we stress that such compar- 
isons can only be made, provided that other observables 
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FIG. 28: (Color on-line) The single-particle source sizes in- 
ferred from the non-identical simulations versus centrality. 
Closed points are for LCMS sizes, open (shifted for clarity) 
for PRF. Circles are pious, squares: kaons, stars: protons. 
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FIG. 29: (Color on-line) Emission asymmetry consistency 
cross-check versus centrality. Closed points are for all par- 
ticles, open (shifted for clarity) for primordial only. 



related to source dynamics, such as e.g. particle spectra, 
elliptic flow and identical particle femtoscopy are repro- 
duced in the model. These tests have been performed for 
Therminator in earlier works [1^ H^] . 

We have calculated the single-particle source sizes for 
all centralities and plotted them in Fig. [551 An ordering 
of source sizes with particle mass is seen as expected, 
as well as increase of the overall source size with Npart- 
Sizes for PRF are larger, as the out radius is scaled by 
the (jt) factor. 

As a final cross-check we show in Fig.l^Hthat the emis- 
sion asymmetries do add up to zero within the systematic 
error of the measurement. 



VIII. CONCLUSIONS 

The THERMINATOR+Lhyquid model was briefly intro- 
duced, and shown to have a unique set of features: hydro- 
dynamic expansion and inclusion of resonance propaga- 
tion and decay for all known particle types. Both features 
are important for non-identical particle femtoscopy. Sim- 
ulations with this model provided predictions for emis- 
sion asymmetries between particles of different masses. 



and validated the analysis methods of non-identical par- 
ticle femtoscopy. 

We underlined the importance of the emission asym- 
metry measurement and provided predictions for RHIC 
energies. We showed that "emission asymmetry", or the 
non-zero difference in the mean emission points of par- 
ticles of same velocity but different mass is predicted to 
arise in hydrodynamical calculation. It was shown to 
be intimately related to the collective behavior of matter 
and to be a direct consequence of the x—p correlations. It 
is the most direct and unambiguous femtoscopic signal of 
such behavior, and provides an independent strong con- 
straint on the models aiming to describe the space-time 
evolution of the heavy-ion collision. Specific predictions, 
both qualitative and quantitative, have been given for 
such asymmetries, for all the considered pair types (pion- 
kaon, pion-proton, kaon-proton) . It was shown that the 
lighter particle is always emitted closer to the center of 
the source, giving a negative emission asymmetry fiout in 
the out direction. Predictions for the size of the source, 
including the dependence on centrality, were also given, 
which were cross-checked with the more precise results 
from identical particle femtoscopy. 

The influence of strongly decaying resonances on emis- 
sion asymmetries was also studied in detail. It was shown 
that due to a specific combination of decay kinematics' 
properties for particles of different masses, the emission 
asymmetry produced by an earlier hydrodynamic stage is 
further magnified. At the same time the resonance decay 
process alone produces emission asymmetry an order of 
magnitude smaller than the x—p correlation, so it cannot 
be used as an alternative explanation of the phenomenon. 

We have presented, in some detail, the theoretical ba- 
sis and some technical aspects of the non-identical parti- 
cle femtoscopy, to be used in ultra-relativistic heavy-ion 
collisions. In particular it was shown how to access the 
"emission asymmetry" , a piece of femtoscopic informa- 
tion that can only be accessed via this type of analysis. 
It was pointed out that spherical harmonics representa- 
tion of the correlation function has useful synergies with 
the analysis; the use of such representation was strongly 
advocated. Important corrections to the correlation func- 
tion were identified: pair purity and momentum resolu- 
tion, and experimentalist recipes were given for apply- 
ing them to the real data. Extracting femtoscopic infor- 
mation from the non-identical particle correlation func- 
tion requires a numerical fitting procedure. It was de- 
scribed in detail, and implemented as a computer code. 
It was also shown that, after application of the proper 
purity correction, the fit procedure was able to recover 
the "true" model input values, which was an important 
validation of the method. 

An estimate of the "fraction of correlated pairs" , com- 
ing from significant non- Gaussian shape and long-range 
tails in the pair separation distributions, has been ob- 
tained, based on the THERMINATOR simulation. The value 
can be an input to an experimental analysis of RHIC 
data, however one must keep in mind that it is model 
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dependent. The systematic error coming from this theo- 
retical estimate was also given. 

Finally, the "experimentalists" analysis procedure, de- 
scribed in the paper, was applied to the correlation func- 
tions obtained from the model calculations. It was shown 
that the method reliably recovers the input model values 
in realistic conditions. Internal consistency cross-checks 
were proposed and tested. The way to compare results 
of non-identical particle femtoscopy and the femtoscopic 
sizes from identical particle analysis was presented; the 
THERMINATOR predictions were shown to be in agreement 
with available identical particle HBT radii results from 
RHIC. 
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Appendix A: Emission function measures 

In femtoscopy one considers pairs of particles. Of par- 
ticular interest are observables related to their relative 
separation r. We refer the reader to the beginning of 
Section mil for the description of symbol conventions, the 
reference system and the relevant reference frames. The 
transition from LCMS to FRF is simply the boost along 
the outwards direction, with the transverse velocity of 
the pair /3t: 

Kut = It {rout - PtAt) 

'^side ~ '^side 
'^long ~ '^long 

At* - jt{At-Ptrout). (Al) 

However, in our calculation we always use the equal time 
approximation, which means that we neglect the time 
difference At* in the PRF. From the components we also 
calculate the length of the relative separation vector: 



side 



long 



(A2) 



The two-particle emission function is, from a mathe- 
matical point of view, any function of the separation four- 
vector, S{Ax), where Ax = {At, r} can be expressed in 
any reference frame. In femtoscopy we use specific func- 
tional forms of S, which are characterized by parameters. 
One must keep in mind that even though the relation be- 
tween relative separation r values in LCMS and PRF is 
given by simple equation Eq. (jAll) , the relations between 
source function parameters (_R's and fi) in these reference 
frames may not be so simple. Moreover, even if there ex- 
ists a simple parametrization of the source in 3D variables 



{rout, fside, Tiong), it does not mean that this source dis- 
tribution expressed in the magnitude of the relative sepa- 
ration r will also have a simple parametric form. The aim 
of this Appendix is to derive relations between commonly 
used source function parameters, defined in various refer- 
ence frames as well as between 3D and ID representation 
parameters. Similar considerations, for different source 
parametrizations have been done in the past [92]. 

We must also repeat the description of Eq. (l2Tt where 
we state that the two-particle emission function Sab is 
(in case of independent emission) a convolution of two 
single particle emission functions. Let us emphasize this 
complication of the femtoscopic measurement: we are in- 
terested in the single particle emission functions, as they 
directly characterize the source. On the other hand fem- 
toscopy provides information about the two-particle (or 
relative) emission functions only. Moreover, the convolu- 
tion procedure of Eq.(|2T|) is not reversible. It means that, 
in a strict mathematical sense, one cannot recover the full 
information about the individual single particle emission 
functions by measuring only the separation distribution. 
Nevertheless, in femtoscopic measurements, reasonable 
additional assumptions can be made and we are able to 
recover some of the information. In this Appendix we will 
discuss how to extract the single-particle emission func- 
tion parameters from the measured two-particle ones. 



1. Traditional emission function parametrizations 

We list here the traditional functional forms of the 
emission function that are used in femtoscopy and list 
their parameters. The easiest is the 1-dimensional ''Rinv" 
parametrization: 
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r exp 
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(A3) 



Note that the formulas neglect the At* dependence, 
which is possible thanks to the equal-time approxima- 
tion in PRF, allowing one to neglect At* dependence of 
I^E'I^. The change from the 3D variables to the ID one 
requires the introduction of the proper Jacobian. Let us 
also explain the factor of 4 before the "Gaussian sigma" 
Rinv instead of the usual 2. It is a particular property of 
a Gaussian distribution that a convolution of two Gaus- 
sians is also a Gaussian with its sigma being the quadratic 
sum of the sigmas of the individual distributions. The 
Rinv parametrization is used for identical particles, for 
which Eq. (pij) can be simplified by noticing that both 
emission functions are the same. If ones assumes that 
a single-particle emission function is a Gaussian, then 
the two-particle one is also a Gaussian with a multi- 
plied by ^/2. Therefore, by fitting the two-particle dis- 
tribution with the functional form (jA3|) one conveniently 
obtains the single-particle a: the Rinv Therefore, the 
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Rinv "HBT radii" reported by experiments are the single- 
particle emission function Gaussian widths, where the 
emission function is assumed to be a 3D Gaussian in the 
PRF. 

A more sophisticated emission function form, used by 
all RHIC and SPS experiments to report identical pion 
femtoscopy results is: 



Sadir) = exp 



'^out 



^ side 



(A4) 



We again note the factor of 4 before the cr's, making them 
the single-particle widths. The important difference is 
that this emission function has three independent width 
i?o„t, and Rfong they are defined in the LCMS 

(hence the L superscript), not in PRF. Note that, un- 
like in PRF, we do not use equal-time approximation in 
LCMS. For identical particles this is not a problem: only 
the Coulomb part of the wave-function depends on At, 
but it is factorized out in the usual fitting procedures, 
and replaced by the averaged value. For non-identical 
particles one can adopt one of two approaches. One can 
continue to use Eq. (jASP to be consistent with identi- 
cal particle femtoscopy, but in this case the Ai spread 
will get absorbed into the three spatial radii. Or one 
can add the fourth component with the time spread to 
Eq. (jASp . Since we aim to compare the results from 
non-identical and identical particle analysis and since we 
are in any case only able to recover the dimensionally- 
averaged source size for non-identical particles we chose 
the first solution. 

For non-identical particles we must make one addi- 
tional modification, the mean emission point should be 
allowed to differ from zero, at least in the out direction. 
Moreover the simple connection between the single- and 
two-particle sizes is no longer possible, since the under- 
lying single-particle emission functions are now different: 



53^(7-) = exp (- 
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(A5) 

We introduced an additional parameter: the mean of the 
distribution in the out direction /i^^^ and the source sizes 
are now the two-particle ones. 

In order to facilitate the comparison between ID and 
3D source sizes in LCMS we also introduce the emission 
function with one, directionally averaged, source size in 
LCMS RL: 



exp 



(A6) 



2. Relating ID and 3D source sizes 



The non-identical particle femtoscopy has been, so far, 
limited to the measurement of the directionally averaged 



source size Rav A question arises, what is the relation 
between R^^ and three-dimensional source sizes i?o«t, 
Rside and Riong- Let us write explicitly the form of the 
emission function in magnitude of r: 
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To find the ID source size corresponding to the 3 3D 
ones, we assume that the above distribution can be ap- 
proximated by: 



r exp 



(A8) 



First, let us note that Sf^^ simplifies exactly to S^^ only 
in the special case of R^ut = R'^ide = Rfong^ which case 
R^y = Rout- If this condition is not met, the S*!^ is not a 
Gaussian in r, and the exact mathematical relation be- 
tween the ID and 3D sizes does not exist. Nevertheless, 
for realistic values of radii, Si^^ is not very different from 
a Gaussian and can be well approximated by 5";^^. One 
can find an effective approximate relation between R^^ 
and {RoutJ Rside^ Rfong) numerically in the following way. 
One generates a significant sample of triplets (rout, rside, 
fiong) where rout is randomly generated from a Gaussian 
of width Roun '''side with width etc. Then, one con- 
to which 



structs a distribution r^^t + 



+ r 



long I 



one then fits numerically the functional form (|A8I) . The 
results of the fit, is the approximate ID source size 
that we seek. We may also want to restrict the fit range 
to small values of r in order to minimize the dependence 
on non-Gaussian features which will mostly affect the 
large-r region. We have performed such calculations and 
concluded that the approximate relation is: 
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)/3. (A9) 



Note that the pair velocity does not enter into the deriva- 
tion, so it is equally valid for LCMS and PRF. 

In addition, in non-identical particle femtoscopy one 
is able to access the first moments of the source 
distribution. Then, S has the general form of: 

S{Rout^ Rside^ Rtng^ l^out, l^side^ l^tng)- ^1^^ tO Symme- 
try relations /if'j^jg must vanish. For collider experiments 
with symmetric rapidity acceptance /i^„g vanishes as 
well, which leaves one additional parameter ii^uf By 
performing a numerical procedure very similar to the one 
in the previous paragraph (the only difference being that 
the rout is now randomly generated from a Gaussian with 
the mean of Hout) obtained an equivalent approximate 
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effective relation: 



R„ 



Ri 



1/5- 



0.3/.i, 



1.3(7/. With these values we write the set of equations 
for the particle source sizes: 



(AlO) 



3. Relating single-particle and two-particle sizes 

As already mentioned for identical particles the rela- 
tion between two-particle and single particle sizes is a 
trivial factor if the assumption of a Gaussian emis- 
sion function is made. For non-identical particle sizes 
such simple connection is not possible, and even an ap- 
proximate one can only be made after certain simplifica- 
tions are done. Let us assume that the two particle types, 
A and B are emitted according to the emission functions 
Sa and Sb which are Gaussians. The two particle emis- 
sion function is then: 



SABir) 



exp 



X exp 



Sirout 

Sir side 



^ ( {■'•'out f^out. 


? 


A 2 
•^side 


^A 






T^A 2 


9 


- uB )2 

\-^out r'out) 




B 2 
■^side 


•^long \ 


2RB ^ 






2R^ ^ / 

■^^long / 


_ r^A , B \, A 1 B 





long 



■^long 



)'^-^long 



ide ""^side 
"'•^long- 



Performing the integration and neglecting the unimpor- 
tant normalization constants one obtains: 



Sab = exp 



X exp 



X exp 



'^{-^out + ^out ) 
side I 



long 



.'^i^tong^ + Rfong^) 



(A12) 



giving us immediately ^^^^ — /i^ 



f^out 



and R^^ = 



R^ + R^ , as expected. Obviously, one cannot re- 
cover the two single-particle source sizes from the one 
two-particle size. However, in this work we have cal- 
culated three independent two-particle sizes: for pion- 
kaon, pion-proton and kaon-proton systems. Therefore 
we have a set of three equations (|Al2p . with A and B 
being tt and K in the first one, tt and p in the second one 
and K and p in the third one respectively. This set of 
equations has three unknowns, the single particle source 
sizes: i?", i?^^ and Rp, so we can solve it to calculate 
them. 

The procedure has to be carried out in a few steps. We 
start with the fit values a/. From this we calculate the 
approximate overall averaged source size -R™, following 
our fit assumptions: Rout = f/i Rside = E^nd Riong = 



„TrK 

a f = 



Kp 



af + af 



(A13) 



and solve it obtaining: 



-J = 












7rp2 

a/ + 








+ a/ 





(A14) 



From these fit parameters one then calculates direction- 
ally averaged single-particle radius R^^ with Eq. (jA9l) - 
cr/ is treated as R^ut- These can be compared to the 
directionally averaged radii in LCMS, which can be cal- 
culated with Eq. (|A9[) from the single-particle 3D HBT 
radii published by the RHIC experiments. 



(All) 



4. Dealing with non-Gaussian source functions 

Let us now consider a non-Gaussian source, as it is ev- 
ident, e.g. from this work, that realistic models predict 
that the emitting source will not have a perfect Gaus- 
sian shape. However, the two-particle emission function, 
defined as a convolution of two single-particle sources, is 
required to be closer to Gaussian than the single-particle 
functions. And the single particle functions themselves, 
even though non-Gaussian, are not very different: they 
usually have a large peak and some long-range tails. 
From this we conclude that even though we know that 
the two-particle source function is not Gaussian, we may 
assume that it will have Gaussian-like features: it will 
have a peak, and the distribution around that peak can 
be reasonably approximated by a Gaussian. It will also 
usually have long-range non-Gaussian tails. This means 
that we can define two variables that can be used as mea- 
sures of the variance: either the RMS of the distribution, 
or the cr of a Gaussian fit to the peak structure in the 
emission function. For the study of non-identical particle 
averaged sizes a is the proper variable to be used. That is 
because pair wave functions, which are the source of fem- 
toscopic correlation, contain strong structures near zero 
r* , while at large r* they do not produce any correlation. 
On the other hand the RMS of the distribution will be 
very sensitive to the long-range tails, which is undesir- 
able in our case. Therefore we define that a "femtoscopic 
size" R for a non-Gaussian emission function is simply 
the (7 of a Gaussian fitted to the two-particle emission 
function near its peak. 
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5. Obtaining Rav from models 



In models we know the emission function exactly. 
Therefore we can infer Rav directly from it. We do this 
numerically, in a manner very similar to the one shown 
in previous sections. First we determine the probabil- 
ity distributions for various components of r: fo{fout)i 
fs{rside), fiiriong)- For the THERMINATOR Calculation, dis- 
cussed in this work, we have found that fo has the form 
of (pS]) , fs is Gaussian, and /; has the form (|37)) . Now we 
proceed with the generation of triplets, construction of 
the S{r) and numerical fitting to obtain Rav Such Rav 
can then be directly compared to the Rav obtained from 
the "experimental" fitting procedure described above. 



R*out = Rout {it) 

R-side ~ Rside 

R-long ~ Rlong 

IJ-out = {it) IJ-ouf 



(A15) 



6. Relations between source sizes in LCMS and 

PR-F Using Eqs. (1X9)) . (jAlOp one can obtain the relations 

for the averaged ID radius as well. One must also 
We can write approximate relations between source independently determine the (74) for the pair sample used 
sizes in the two reference systems, they are: in a given analysis. 
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